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NEAR  FIELD  SOURCE  CHARACTERIZATION 
OF 

CONTAINED  NUCLEAR  EXPLOSIONS  IN  TUFF 


BRIAN  W.  STUMP 

Air  Force  Weapons  Laboratory 
Klrtland  AFB  NM  87117 


ABSTRACT 


The  seismic  source  of  three  contained  nuclear  explosions  In  tuff  at  Pahute 
Mesa  is  characterized  by  applying  the  moment  tensor  source  representation  to 
near-field  accelerometer  recordings  of  the  events.  Propagation  path  effects 
are  accounted  for  In  the  study  by  generalized  ray  calculations  for  several 
plane  layered  models.  Synthetic  seismograms  generated  from  the  predicted 
sources  model  well  the  distance  and  azimuthal  variations  In  the  data.  The 
moment  tensors  are  dominated  by  the  Isotropic  component  with  the  absolute 
ratio  of  isotropic  to  devlatorlc  component  varying  between  2  and  20,  depending 
on  the  partlcualr  Green's  function.  The  Isotropic  sources  Indicate  the  proce¬ 
dure  Is  a  good  measure  of  relative  yield  while  the  absolute  yield  Is  dependent 
upon  the  exact  Green's  function.  Source  time  functions  Indicate  the  possibi¬ 
lity  of  a  double  pulse  source  separated  by  1.0  to  1.5  seconds.  This  Interpre¬ 
tation  Is  bounded  by  the  fact  that  phase  Information  Is  the  least  well  resolved 
In  the  study.  In  contrast,  the  stability  of  the  spectral  characteristics, 
such  as  corner  frequency,  high  frequency  roll-off,  and  peak  value  suggest  the 
use  of  the  moment  tensor  In  the  frequency  domain  for  yield  and  discrimination 
studies. 


SOURCE  CHARACTERIZATION 


Given  a  set  of  observed  waveforms  from  an  explosion  or.  an  earthquake  and 
knowing  the  propagation  path  effects  of  the  geological  material,  a  quan¬ 
titative  procedure  for  determining  the  source  function  In  space  and  time  Is 
desirable.  The  moment  tensor  source  representation  Is  such  a  formulation 
(Gilbert,  1970;  Gilbert  and  Dzlewonskl,  1975;  Backus  and  Mulcahy,  1976  a  and 
b;  Stump  and  Johnson,  1977;  Backus  1977  a  and  b;  Strelltz,  1977). 

Assuming  an  explosion  or  earthquake  source  can  be  represented  as  a  set  of 
equivalent  body  forces,  then  the  source  can  be  written  as  a  series  of  moments. 
For  small  sources  or  large  wavelengths,  only  the  first  term  of  the  series  Is 
retained,  and  the  displacement  at  any  point  and  time  can  be  written  as: 


Uv'i'.f)  =  Gw,;.  u‘,+.\0,o) 


Mylfl.-t') 


Where  (1^  Is  the  displacement  In  the  k  direction,  G^  Is  the  Green's  function, 
Mjj  Is  the  moment  tensor,  j  Indicates  derivative  with  respect  to  xj  and  ® 
represents  temporal  convolution.  A  more  complete  derivation  of  (1)  Is  given 


In  Stump  and  Johnson  (1977). 


In  the  frequency  domain,  the  equation  reduces  to: 


,o,o)  .  (2) 

Knowing  the  propagation  path  effects  (6^ t j) ,  one  can  determine  the  source 
(Hjj)  from  a  set  of  observational  data  (Uj<). 

In  this  study,  the  moment  tensor  formulation  will  be  applied  to  three  sets 
of  data  from  contained  nuclear  explosions.  The  sum  of  the  diagonal  elements 
of  the  moment  tensor,  the  Isotropic  component,  is  proportional  to  volume 


changes  in  the  source  region  and  will  be  Interpreted  as  the  explosion  com¬ 
ponent  of  the  source.  This  component  can  then  be  compared  to  the  remaining 
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portion  of  the  moment  tensor,  Its  devlatorlc  component,  to  analyze  that  por¬ 
tion  of  the  source  not  accounted  for  by  a  symmetric  explosion.  * 

Mtr  (Isotropic)  ■  1/3  M-f  j  Si j  (3) 

Dij  (devlatorlc)  -  Mij  -  Mjr  Si j  (4) 

If  one  has  some  Idea  of  the  source  characterization,  then  constraints  can 
be  placed  on  Equation  2.  One  such  constraint  Is  that  the  source  Is  a  sym¬ 
metric  explosion,  but  the  time  function  Is  unknown.  Equation  2  reduces  to: 

Uk(x*.f)  *  M0(Sklfl  +  6|t2>2  +  Sk3,3)  F(f)  <5) 

F(f)  Is  the  unknown  source  time  function  In  the  frequency  domain. 
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DATA  ANALYSIS 

This  report  will  discuss  the  analysis  of  near-field  accelerometer  data 
from  three  nuclear  explosions  detonated  In  Pahute  Mesa  at  the  Nevada  Test 
Site.  The  three  shots  studied  are:  Handley  with  a  Wood-Anderson  magnitude 
determined  by  the  U.C.  Berkeley  seismographic  station  of  6.3,  Pipkin  with  a 
Wood-Anderson  magnitude  of  5.5,  and  Farm  which  is  similar  in  size  to  Pipkin. 

The  instrumentation  arrays  for  Handley  and  Pipkin  are  illustrated  in 
Figure  1.  Eight  3-component  accelerometers  were  fielded  over  a  180°  arc  at  8 
kilometers  from  the  Handley  event.  Stations  1,  2,  4,  and  5  operated  properly 
yielding  twelve  components  of  motion.  A  similar  array  to  Handley  was 
installed  for  the  Jorum  event.  While  the  equipment  was  operating,  the  Pipkin 
explosion  was  detonated  yielding  data  at  a  variety  of  azimuths  (80°  arc)  and 
ranges  (2-14  km).  Stations  1,  2,  3,  4,  5,  and  6  operated  for  a  total  of  14 
components  of  motion  (not  all  accelerometers  functioned  properly). 

The  Farm  Instrumental  array  Is  given  in  Figure  2.  The  3-component  acce¬ 
lerometers  spanned  135°  In  azimuth  and  varied  in  range  between  2  and  10 
kilometers.  Of  the  six  stations,  only  four  operated— 2,  4,  5,  and  6— yielding 
twelve  components  of  data. 

Figure  3  illustrates  the  three  Instrumental  arrays  with  respect  to  one 
another  and  the  Silent  Canyon  Caldera,  the  major  structural  feature  of  Pahute 
Mesa.  The  Pipkin  and  Farm  events  are  both  within  the  caldera,  while  Handley  Is 
out.  The  accelerometers  for  Farm  are  all  within  the  caldera,  while  those  for 
Pipkin  and  Handley  Farm  are  inside  and  outside. 

The  Instruments  used  in  all  three  of  these  experiments  were  force-balance 
accelerometers  with  a  flat  response  between  0.02  to  50  Hz.  The  outputs  of  the 
gages  are  5  volts/g  for  the  verticals  and  7.5  volts/g  for  the  horizontals. 
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For  Pipkin  and  Handley,  the  data  was  recorded  on  analog  tape  at  NTS.  The 
seismograms  were  low  pass  filtered  at  10  Hz  with  a  4-pole  Butterworth  filter. 
The  data  was  then  digitized  at  a  rate  of  54  samples/secondl 

Digital  recording  of  the  data  for  the  Farm  test  utilizing  Sprengnether 
DR-100  digital  event  recorders  was  accomplished.  The  data  was  filtered  with  a 
2-pole  low  pass  Butterworth  at  30  Hz  and  a  2-pole  high  pass  Butterworth  at  5 
seconds  prior  to  digitization  at  100  samples  per.  second. 

Analysis  of  the  data  and  noise  samples  in  the  frequency  domain  indicates 
that  there  are  problems  with  the  long  period  data  between  2  and  5  seconds  from 
all  three  events.  This  noise  becomes  increasingly  evident  as  one  tries  to 
calculate  velocities  and  displacements.  Since  the  source  studies  using  this 
data  will  essentially  deconvolve  the  Instrument  and  propagation  path  effects 
out  of  the  data,  one  must  be  careful  with  the  data  at  long  periods.  In  order 
to  avoid  long  period  problems,  some  post  processing  of  the  observational  data 
was  done.  First  a  high  pass  Butterworth  filter  was  applied  to  the  data.  The 
slope  and  corner  frequency  of  the  filter  was  determined  from  the  noise 
estimate,  a  desire  to  eliminate  the  rise  of  the  spectra  at  long  periods,  and 
an  attempt  to  yield  realistic  velocity  and  displacement  records  free  of  ramps, 
steps,  and  very  long  period  noncausal  signals.  Typically,  the  filter  corners 
were  between  2  and  5  seconds  with  1  or  2  slopes  at  high  frequencies.  After 
application  of  these  filters,  the  data  were  corrected  for  a  DC  offset  in  velo¬ 
city  or  a  ramp  In  displacement. 

It  Is  this  corrected  data  In  velocity  and  acceleration  which  were  used 
for  the  source  study.  The  corrected  acceleration,  velocity,  and  displacement 
records  from  the  Handley  event  are  given  In  Figures  4  and  5.  The  uncorrected 
displacement  spectra  for  each  observation  are  given  In  Figures  6  through  17. 
The  rise  of  the  spectra  at  long  periods  can  be  noted  In  these  figures.  The 


4 


time  domain  data  yields  vertical  components  of  motion  which  are  quite  correc¬ 
table  while  the  coherence  of  the  radial  and  transverse  motions  Is  less.  The 
transverse  motions  in  acceleration  are  comparable  in  size  to  the  vertical  and 
radial  In  many  cases.  In  velocity  and  displacement,  the  transverse  motions 
are  on  the  average  comparable  in  size  to  the  radial  motion  and  half  the  size 
of  the  vertical. 

Unlike  the  Handley  data  which  were  all  at  the. same  radial  distance,  the 
Pipkin  data  span  distance  2  to  14  km  from  the  source  (Figure  1).  The  data 
were  post  processed  in  the  same  manner  as  that  for  Handley.  The  corrected 
accelerations,  velocities,  and  displacements  for  Pipkin  are  given  In  Figures 
18  through  20.  The  displacement  spectra  prior  to  correction  are  given  in 
Figures  21  through  35.  The  frequency  content  of  the  waveforms  apparently 
Increases  with  range.  This  observation  Is  quite  apparent  when  comparing  the 
vertical  acceleration  at  Station  1  (2.77  km)  and  the  vertical  acceleration  at 
Station  6  (13.55  km)  (Figures  18  and  20).  Again,  the  relative  ratios  of  the 
radial,  transverse,  and  vertical  motions  are  similar  to  those  observed  for 
Handley. 

The  Farm  data  was  the  only  digitally  recorded  set.  Again,  the  data  were 
post  processed  In  the  same  manner  as  the  previous  two  data  sets.  The 
corrected  ground  motion  data  are  given  in  Figures  36  and  37  with  the  accom¬ 
panying  spectra  In  Figures  38  through  49.  The  data  set  has  transverse  motions 
which  are  comparable  to  the  radial  and  vertical.  The  accelerograms  Indicate  a 
good  deal  of  complexity.  There  Is  the  Indication  of  a  fairly  large,  long 
period  secondary  arrival  on  the  displacement  records  at  approximately  four 
seconds.  This  arrival  is  particularly  apparent  on  the  vertical  components  of 
Stations  4  and  6  (Figure  36). 

Summarizing  the  observational  data,  we  have  good  quality  accelerometer 
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data  from  three  explosions  that  have  been  corrected  for  long  period  noise. 

The  accelerations  yield  transverse  motions  which  are  comparable  to  the  radial 
and  vertical  records.  At  longer  periods,  as  exemplified  by  the* velocity  and 
displacement  records,  the  vertical  motions  are  on  the  average  a  factor  of  two 
larger  than  the  radial  and  transverse.  The  data  appear  to  increase  in  fre¬ 
quency  with  range  over  2-14  km.  For  data  at  the  same  range,  the  vertical 
records  from  varying  azimuths  are  quite  correctable,  while  the  radial  and 
transverse  motions  are  much  less  so. 


PROPAGATION  PATH  EFFECTS 


Utilizing  geological  maps,  geophysical  explorations*  and  the  observational 
ground  motions  from  the  three  explosions*  a  propagation  path  model  for  Pahute 
Mesa  at  NTS  was  developed.  Even  though  there  are  great  geological  variations 
across  Pahute  Mesa— water  table  fluctuations,  terrain  differences,  and  struc¬ 
tural  boundaries,  such  as  the  Silent  Canyon  Caldera  (Figure  3)— It  was  our 
desire  to  check  the  feasibility  of  producing  an  average  plane  layered  model 
for  the  mesa. 

The  construction  of  the  model  began  with  the  early  work  of  Hamilton  and 
Healy,  1969.  Their  Initial  velocity  model  for  the  Benham  explosion  on  Pahute 
Mesa  was  supported  by  a  4.1  km  deep  well  log,  geologic  data,  and  long  line 
refraction  studies  of  NTS.  Their  model  Is  reproduced  in  Table  1.  Although 
this  model  Is  good  at  deep  depths,  we  were  interested  in  the  fine  structure 
above  one  km  since  two  of  the  shots  studied  were  emplaced  In  this  depth  range. 

Figure  50  is  the  well  log  above  the  Farm  working  point.  The  acoustic  log 
Indicates  a  gradual  increase  In  velocity  with  depth  from  about  2.0  to  3.4 
km/sec  with  the  possibility  of  a  velocity  jump  at  the  water  table.  One  must 
keep  in  mind  that  the  logging  method  changes  at  the  water  table.  Using  this 
data  and  that  from  additional  drill  holes  at  NTS,  an  initial  velocity  model 
was  developed.  The  model  Is  very  close  to  that  of  Helmberger  and  Hadley 
(1981). 

The  complexity  in  the  velocity  model  of  Interest  is  controlled  by  the 
wavelengths  of  energy  leaving  the  explosive  source.  Although  the  velocity  log 
of  Figure  50  shows  great  variability  on  a  local  scale,  we  wish  some  type  of 
average  structure.  In  order  to  help  further  constrain  the  Pahute  Mesa  velo¬ 
city  model  for  the  wavelengths  of  the  nuclear  sources,  the  accelerometer  data 
from  the  Farm  and  Pipkin  experiments  were  used  to  obtain  arrival  time  data. 


TABLE  1 


i  to  Laver  (km) 

Lav  -  Velocity  (km/sec). 

0.0 

2.7 

0.94 

3.4 

1.33 

3.8 

2.14 

4.4 

2.50 

5.1 

5.00 

6.1 

25.00 

7.0 

35.00 

8.0 
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This  time  data.  In  conjunction  with  travel  time  curves  from  the  Pahute  Mesa 
model,  were  utilized  to  refine  the  structure.  Farm  and  Pipkin  were  detonated 
at  similar  depths,  but  Pipkin  was  just  above  the  water  table  whfle  Farm  was 
just  below  it.  Since  the  water  table  Is  associated  with  a  velocity  jump  In 
our  model,  the  travel  time  curves  for  the  Pipkin  and  Farm  shots  are  quite  dif¬ 
ferent,  and  the  effect  must  be  taken  into  account  when  analyzing  the  travel 
time  data.  The  finalized  velocity  model  structure,  along  with  the  explosion 
locations,  is  given  In  Figure  51.  The  travel  time  curves  for  a  source  just 
above  the  water  table  with  the  Pipkin  arrival  time  data  are  given  in  Figure 
52.  The  travel  time  curve  for  a  source  below  the  water  table  with  the  Farm 
arrival  time  data  Is  given  In  Figure  53.  The  observational  data  fit  the  pre¬ 
dictions  of  the  models  quite  well.  The  change  in  crossover  distance  from  5 
km  to  9.5  km  as  one  moves  down  across  the  water  table  at  0.65  km  seems  to 
match  the  observed  arrival  time  differences  for  the  two  shots.  In  the  ray 
tracing,  the  second  interface  at  1.45  km  was  Ignored  as  we  feel  the  layering 
at  this  point  replicates  a  gradient  In  velocity.  Inclusion  of  this  Interface 
as  a  refraction  led  to  predicted  arrivals  before  those  observed  In  the  actual 
arrival  time  data. 

The  Interpretation  of  the  proposed  structure  Is  that  the  0.65  km  Interface 
represents  the  water  table  depth  and  Is  constrained  by  well  levels  In  the 
caldera.  As  one  moves  out  of  the  caldera  near  the  Handley  site,  the  water 
moves  to  half  this  depth.  Since  we  are  Interested  in  an  average  velocity 
model  and  that  the  ray  paths  to  the  Handley  gages  propagate  both  In  and  out  of 
the  caldera,  the  average  depth  of  0.65  km  to  the  water  table  will  be  used.  As 
we  already  discussed,  the  Interface  at  1.45  km  replicates  the  velocity  gra¬ 
dient  with  a  two-layer  structure.  Finally,  the  2.50  km  Interface  represents 
the  turning  point  of  the  downgoing  energy  In  the  structure.  Its  depth  and 
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velocity  are  constrained  by  the  arrival  time  data  beyond  6  km.  We  feel  the 
essential  elements  of  this  model  are  the  water  table,  the  velocity  gradient 
below  the  source,  and  the  turning  point  of  the  rays  at  2.50  km. 

Up  to  this  point,  the  focus  has  been  upon  the  compresslonal  velocity 
structure.  The  constraints  upon  the  shear  velocity  are  much  less.  As  we 
shall  see  later,  the  shear  structure  Is  Important  In  the  Green's  functions  for 
they  control  the  ratio  of  the  energy  from  Isotropic  to  devlatorlc  sources. 

The  shear  wave  velocity  also  controls  the  Importance  of  P  to  S  conversions  at 
the  Interfaces.  The  shear  velocities  were  determined  from  both  arrival  time 
data  and  frequency-wave  number  analysis  of  small  array  data  recording  explo¬ 
sion  data  at  Pahute  Mesa  (K.  McLaughlin,  U.C.  Berkeley,  personal  communication). 
TheyS-news  model  (Figure  51)  Involved  material  with  an  average  Poisson's  ratio 
of  0.28.  Since  we  are  very  interested  In  the  ratio  of  P  to  S,  a  slow  shear 
model-y3-olds— was  also  Investigated.  This  model  began  with  a  Poisson's 
ratio  of  0.45  In  the  surface  layer  and  slowly  decreasing  to  0.26  In  the 
half-space. 

Using  the  velocity  structure  of  Figure  51,  synthetic  Green's  functions 
were  calculated.  Models  utilizing  both/3-news  and/S-olds  were  Included. 
Synthetics  were  run.  Including  turning  rays  and  conversions  at  Interface  two 
and  not  Including  them.  The  Green's  functions  were  run  for  a  source  just 
above  and  just  below  the  water  table  to  replicate  the  Pipkin  and  Farm 
emplacements.  All  the  synthetics  for  the  layered  structure  were  calculated 
using  the  generalized  ray  method  (Helmberger,  1968;  Pao  and  Gajewskl,  1977). 

For  comparison,  half-space  Green's  functions  were  also  calculated  (Johnson, 
1974).  The  half-space  propagation  path  was  Included  because:  (1)  It  Is  the 
simplest  of  all  path  corrections,  (2)  the  results  Include  near-field  effects, 
and  (3)  Rayleigh  wave  contributions  are  Included. 


The  Green's  functions  for  a  source  just  below  the  water  table  (Farm, 

Figure  51)  are  given  In  Figures  54  through  65.  We  will  use  this  source  con- 

« 

figuration  to  investigate  the  various  velocity  models.  All  the  results  are 
delta  function  responses  passed  through  the  accelerometer.  As  we  noted  In  the 
theory  section,  to  completely  characterize  the  material  response,  ten  sources 
are  needed.  Primary  emphasis  shall  be  placed  on  the  radial  and  vertical 
response  to  an  explosive  source. 

Figures  54  and  55  are  the  radial  and  vertical  explosive  Green's  functions 
for  the  structure  of  Figure  51  ^6-news)  Ignoring  conversions  and  turning  rays 
at  the  second  interface.  Notions  between  2  and  12  kilometers  are  given.  A 
total  of  31  rays  were  included  In  the  calculation.  The  travel  time  curves 
Indicate  a  crossover  distance  between  9  and  10  km  (Figure  S3)  and  this  can  be 
seen  in  the  synthetics.  The  radial  and  vertical  components  are  similar  In 
both  amplitude  and  wave  shape.  The  complete  set  of  ten  Green's  functions  and 
their  spectra  for  a  6.34  km  range  are  given  in  Figure  56.  Radial  and  vertical 
Green's  functions  result  from  all  sources  while  the  transverse  contribution 
comes  only  from  the  strike-slip  and  dip-slip  sources.  Comparison  of  peak  time 
amplitudes  and  spectral  estimates  Indicates  transverse  motions  which  are  13.6 
to  7.6  times  larger  than  the  energy  from  the  explosive  source.  In  order  to 
explain  observed  transverse  motions  equal  In  size  to  the  vertical,  one  needs 
a  source  only  7.6  to  13.6  times  smaller  than  the  explosion.  As  we  will  see, 
this  ratio  varies  for  the  set  of  Green's  functions  tred  In  this  study.  The 
vertical  to  radial  ratio  for  the  explosive  source  for  this  propagation  path 
model  Is  unity. 

In  order  to  Investigate  the  effect  of  the  second  Interface  as  a  turning 
point,  synthetics  were  run  for  the  model  Just  discussed,  but  conversions  and 
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turning  rays  at  the  second  interface  were  Included.  A  total  of  67  rays  were 
calculated.  The  explosion  Green's  functions  at  the  4.20,  6.34,  and  9.67  km 

r 

ranges  are  illustrated  In  Figures  57  and  58.  Comparison  of  Figures  54  and  57 
indicates  that  the  energy  converted  and  turned  at  the  second  Interface  has  a 
significant  contribution  to  the  synthetics.  At  the  9.67  km  range,  the  energy 
off  the  second  Interface  is  actually  the  first  arrival.  The  complete  Green's 
functions  for  the  6.34  km  range  are  given  in  Figure  59.  The  peak  transverse 
motions  are  13.2  to  5.3  times  larger  than  the  explosive  source  motions.  The 
vertical  to  radial  ratio  for  the  explosion  Green's  function  Is  0.8. 

The  effect  of  varying  the  shear  wave  velocity  is  illustrated  in  Figures  60 
through  62.  The  source  structure  is  the  same  as  before,  but  with  the^d-olds 
shear  wave  structure,  larger  Poisson's  ratio  (0.45)  near  the  surface.  The 
radial  accelerations  from  the  explosion  Green's  function  show  greatly 
Increased  P  to  S  conversions,  while  the  vertical  accelerogram  shows  less 
change.  These  conversions  are  most  pronounced  at  the  6  km  or  less  range.  The 
slow  shear  velocity  yields  a  vertical  explosive  Green's  function  which  Is  1.7 
times  larger  than  the  radial.  The  complete  Green's  functions  (6.34  km  range) 
In  both  time  and  frequency  are  given  in  Figure  62.  The  ratio  of  the  trans¬ 
verse  body  waves  to  the  explosion  body  waves  varies  between  13  and  46. 

For  final  comparison,  the  Green's  functions  In  an  elastic  half-space  are 
given  In  Figures  63-65.  The  direct  P  wave  dominates  the  explosive  source 
(Figures  63  and  64).  Only  the  upgolng  energy  Is  returned  to  the  free  surface. 
This  phenomena  Is  unlike  that  of  the  previous  models.  The  layering  In  the 
three  previous  models  could  always  return  energy  to  the  free  surface.  The 
complete  Green's  functions  at  the  6.34  km  range  are  given  In  Figure  65.  The 
seismograms  are  much  simpler  than  before,  consisting  of  P,  SP  diffracted,  S, 
and  Rayleigh  arrivals.  The  vertical  to  radial  ratio  of  the  P  waves  from  the 


explosion  Is  0.36.  The  ratio  of  the  transverse  waves  to  the  explosive  energy 
is  between  1  and  20. 

This  complete  range  of  Green's  functions  will  be  used  In  the  source  Inver¬ 
sion  for  the  Farm  event.  The  effect  of  ranging  from  the  simplest  propagation 
path  model  (elastic  half-space)  to  a  fairly  complex  model  (67-News)  will  be 
shown.  In  particular,  we  wish  to  Investigate  the  sensitivity  of  the  inver¬ 
sions  In  determining  relative  or  absolute  source  strength  and  the  variability 
of  the  Isotroplc/devlatorlc  source  ratios.  The  high  value  of  the  maximum 
transverse  Green's  function  to  the  explosive  Green's  functions  seems  to  Indi¬ 
cate  that  a  primarily  Isotropic  source  may  help  explain  the  large  transverse 
motions  observed  from  contained  explosions. 

For  the  Pipkin  and  Handley  studies,  a  similar  set  of  Green's  functions 
will  be  utilized.  The  Pipkin  source  Is  just  above  the  water  table,  while 
Handley  Is  burled  at  1.2  km,  some  0.55  km  below  our  water  table.  Figure  66 
gives  a  comparison  for  the  explosive  Green's  functions  at  the  Farm,  Pipkin, 
and  Handley  burial  depths  and  observed  at  9.67  km. 
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MOMENT  TENSOR  INVERSIONS 


Utilizing  the  data  and  Green's  functions  already  reviewed  In  Equation  2, 
the  moment  tensor  representations  of  the  Farm,  Handley,  and  Pipkin  explosions 
were  determined.  The  Inversions  were  done  In  the  frequency  domain  and  then 
transformed  back  Into  the  time  domain  to  yield  source  time  functions.  In 
order  to  check  the  adequacy  of  the  resulting  source,  the  moment  tensor,  along 
with  the  appropriate  Green's  function.  Is  substituted  Into  Equation  1  to  yield 
a  calculated  seismogram.  These  predicted  seismograms  are  then  cross  corre- 
lated  with  the  observations  to  give  one  a  simple  measure  of  the  fits. 

In  addition  to  using  the  four  sets  of  Green's  functions  with  each  obser¬ 
vational  data  set,  various  time  windows  were  utilized  Including  the  entire 
waveform  and  just  the  P  arrivals.  Constrained  Inversions  restricting  the 
source  to  purely  Isotropic  were  also  conducted.  A  total  of  six  or  more  Inver¬ 
sions  were  done  on  each  of  the  three  data  sets.  This  large  matrix  of  Inver¬ 
sions  was  done  to  Investigate  the  stability  of  the  procedure  and  Its  sen¬ 
sitivity  to  changes  In  data  and  Green's  functions. 

For  Illustration  purposes,  the  complete  moment  tensor  time  functions  and 
the  synthetic  matches  to  the  observational  data  for  the  Farm,  Handley  and 
Pipkin  explosions  will  be  given  for  only  one  Green’s  function.  The  Green's 
function  that  will  be  focused  upon  Is  the  one  utilizing  the  new  shear  wave 
velocity, yd -news,  and  excluding  conversions  and  turning  rays  at  Interface  two 
(model  32).  The  results  of  the  Inversions  using  the  other  three  Green's  func¬ 
tions  are  sumnarlzed  In  Table  2. 

The  observed  and  calculated  seismograms  for  Farm  are  given  In  Figures  66 
and  67.  A  total  of  10.24  seconds  of  data  was  used  In  the  Inversion.  The 
signal  to  noise  ratio  In  this  data  Is  the  maximum  of  three  sets  and  allowed 
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the  study  of  the  long  data  window.  The  correlation  coefficients  for  the 
observed  and  calculated  seismograms  are  given  In  the  figures.  The  peak  ampli¬ 
tudes  and  the  general  decay  envelope  of  the  data  are  well  modeled.  The  maxi¬ 
mum  correlation  coefficients  are  found  for  the  radial  and  vertical  components, 
0.77  and  0.74,  while  that  for  the  transverse  Is  less,  0.69.  As  a  whole,  these 
fits  are  quite  good. 

The  complete  set  of  near-  and  far-fleld  moment  tensors  from  this  Farm 
Inversion  Is  given  In  Figure  68.  In  addition  to  the  six  components  of  the 
moment  tensor,  the  sum  of  diagonal  elements,  the  Isotropic  moment.  Is  also 
given.  All  moment  tensor  plots  are  scaled  to  the  maximum  of  each  element,  so 
comparison  of  component  sizes  must  be  done  with  the  nunerteal  values  noted  In 
the  figures.  The  moment  tensor  Is  dominated  by  the  three  diagonal  elements  of 
the  tensor  and,  thus.  Its  Isotropic  component  (Equation  3).  These  elements 
all  have  similar  time  functions  and  amplitudes.  The  deviator 1c  components  of 
the  moment  tensor  (Equation  4)  are  a  factor  of  4.4  smaller  than  the  Isotropic 
In  the  far-fleld  and  3.4  In  the  near-field.  The  time  functions  of  the 
devlatorlc  components  do  not  resemble  each  other  or  the  Isotropic  component. 
The  arrow  on  the  Isotropic  time  function  In  Figure  68  Indicates  the  predicted 
source  time  from  the  windowing  of  the  data  and  Green's  functions.  Taking  this 
reference  time,  there  are  two  distinct  parts  to  the  Isotropic  time  function 
separated  by  1.5  seconds.  The  first  part  of  the  time  function  Is  a  single 
pulse  In  the  near-field  with  approximately  2/3  second  duration,  while  the 
second  part  1.5  seconds  later  consists  of  several  cycles  of  0.67  Hz  energy. 
There  Is  little  energy  beyond  5  seconds. 

The  fits  to  the  Handley  data  using  the  same  layered  structure,  but 
accounting  for  the  deeper  source,  are  given  In  Figures  69  and  70.  The  data 
quality  limits  one  to  4.74  second  records.  As  Figure  1  Illustrates,  the  data 
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from  Handley  are  all  at  the  same  range  (8  km).  The  vertical  components  of 

velocity  are  the  most  correctable  between  observations,  and  the  synthetic 

* 

fits  to  this  data  are  the  best.  The  average  correlation  coefficient  for  the 
vertical  data  Is  0.90,  the  radial  data  Is  0.74,  and  the  transverse  data  Is  0.65. 

The  Handley  moment  tensors  are  given  In  Figure  71.  The  arrow  on  the 
isotropic  component  of  the  moment  tensor  Indicates  the  Initiation  time  for  the 
source  as  predicted  by  the  data  and  Green's  function  windows.  The  diagonal 
elements  of  the  moment  tensor  are  all  similar  In  time  and  amplitude,  giving  a 
source  dominated  by  the  Isotropic  moment  tensor.  In  the  far-fleld  time 
functions,  the  Isotropic  moment  tensor  Is  7.67  times  larger  than  the 
devlatorlc  components.  In  the  near-field,  this  ratio  Is  8.33.  The  Isotropic 
source  Is  a  single  pulse  In  this  case  with  a  pulse  width  between  0.67  and  1.00 
seconds. 

The  fits  to  the  final  data  set,  Pipkin,  are  given  In  Figures  72  to  74.  In 
the  past  two  data  sets,  the  vertical  components  of  motion  were  the  most 
correctable.  In  this  data  set,  vertical  components  were  lost  at  two 
stations.  In  addition  to  this  loss,  great  variability  In  frequency  content 
with  range  Is  seen  In  the  data.  A  total  of  4.74  seconds  of  data  are  fit  In 
this  analysis.  The  average  correlation  coefficients  for  the  radial,  vertical 
and  transverse  data  are  0.51,  0.70,  and  0.62  respectively. 

The  ament  tensors  for  this  Inversion  are  given  In  Figure  75.  The  arrow 
on  the  Isotropic  moment  tensor  Indicates  the  theoretical  source  Initiation 
time.  The  amplitudes  of  the  three  diagonal  elements  of  the  moment  tensor  are 
approximately  equal  In  magnitude  but,  unlike  the  previous  two  data  Inversions, 
there  Is  some  variability  In  the  time  functions.  The  ratio  of  the  Isotropic 
to  the  devlatorlc  moment  tensor  In  the  near-field  Is  1.95  and  2.47  In  the 
far-fleld. 
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The  Isotropic  source  spectra  for  the  three  inversions  just  discussed  are 
given  in  Figure  76.  The  solid  lines  are  the  spectral  estimates,  while  the 
dashed  lines  are  variance  estimates.  The  variances  .itersect  the  moment  esti¬ 
mates  between  10  and  20  Hz.  The  maximum  separation  between  the  spectra  and 
the  variances  runs  between  40  to  60  db  near  1  Hz.  The  high  frequency  roll¬ 
offs  for  all  the  sources  are  between  2  and  3.  The  corner  frequencies  vary 
little  between  the  three  sources  with  a  value  between  1  and  2  Hz.  The  spectra 
decrease  at  the  long  periods  yielding  a  peaked  spectra  around  the  corner 
frequency.  Although  not  included  in  this  report,  the  moment  tensor  spectra 
from  inversions  utilizing  the  other  propagation  path  models  yield  spectral 
characteristics  consistent  with  those  reviewed.  Of  all  the  moment  tensor 
characteristics,  their  high  frequency  slopes,  corner  frequencies,  and  maximum 
values  all  remain  Invariant  while  the  phase  information  varies. 

The  isotropic  near-field  moment  tensors  from  the  Farm,  Hand'^y,  and  Pipkin 
Inversions  utilizing  the  three,  layered  Green's  functions  (31,  67,  OLDS)  are 
summarized  in  Figure  77.  The  three  moments  for  Farm  all  Indicate  a  source 
function  with  two  pulses  separated  by  1.0  to  1.5  seconds.  The  period  and 
duration  of  the  second  pulse  decrease  as  one  moves  from  the  F31  to  the  FOLDS 
models.  The  Handley  source  also  Indicates  two  time  functions  separated  1.0 
seconds  In  time.  These  phenomena  are  particularly  apparent  for  the  HOLDS 
model.  There  Is  a  hint  of  two  pulses  In  the  Pipkin  source  too.  The  period  of 
this  source  decreases  as  one  moves  from  the  P31  to  the  POLDS  models. 
Considering  the  tremendous  variations  In  the  observed  data  from  these  three 
explosive  sources,  the  differences  In  propagation  paths,  and  the  physlcal'dlf- 
ferences  In  the  sources  themselves,  the  Isotropic  components  of  the  moment 
tensor  are  remarkably  consistent. 
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The  absolute  size  of  the  moment  tensors  for  all  the  inversions  of  this 
study  are  summarized  In  Table  2.  The  first  observation  Is  that  the  whole 
record  Inversions  (67,  31,  OLDS,  HS)  yield  moments  which  are  nearly  the  same. 
The  Isotropic  moments  change  by  a  maximum  of  52  percent  for  Farm,  28  percent 
for  Handley,  and  121  percent  for  Pipkin.  The  largest  variation  In  the  source 
with  Green's  function  Is  In  the  ratio  of  Isotropic  to  devlatoric  components. 
The  far-fleld  ratio  for  Farm  varies  from  7.33  for  the  F67  Green's  function  to 
27.3  for  the  FOLOS  model.  The  ratio  of  transverse  body  waves  to  explosive 
body  waves  for  the  Green's  functions  varies  from  7.6  to  13.6  for  the  F67  model 
(Figure  56)  to  between  13  and  46  for  the  FOLOS  model  (Figure  62).  As  one 
moves  to  slower  shear  wave  velocities  in  the  propagation  path  models,  the 
Isotropic  to  devlatoric  ratio  In  the  source  function  Increases  due  to 
increased  shear  wave  amplitudes. 

The  moments  of  Farm  and  Pipkin  are  nearly  equal  In  size.  Handley  has  a 
moment  tensor  between  15  and  20  times  larger  than  the  two  other  events. 

Since  the  Initial  P  waves  are  the  most  coherent  arrivals  In  the  data  sets 
and  should  be  least  contaminated  by  shear  waves,  inversions  were  completed 
using  the  first  second  of  the  Farm  data.  This  source  Inversion,  In  conjunc¬ 
tion  with  the  earlier  whole  record  studies,  will  help  us  understand  the  Impor¬ 
tance  of  P  and  S  waves  In  constraining  the  Isotropic  and  devlatoric  components 
of  the  explosive  source. 

The  observed  and  calculated  seismograms  utilizing  the  Farm  data  and  the 
F67  propagation  model  are  given  In  Figures  78  and  79.  The  fits  to  the  radial 
and  vertical  velocities  are  the  best  of  the  trial  inversions.  The  correlation 
coefficients  for  the  radial  velocity  Is  0.82  and  for  the  vertical  velocity  Is 
0.85.  The  transverse  motions  are  poorly  modeled  with  an  average  correlation 
coefficient  of  0.23.  The  Inversion  window  contained  no  shear  arrivals  In 
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Green's  functions,  so  only  near-field  components  of  motion  remained  to  fit  the 
observations.  The  near-field  components  from  this  model  are  not  sufficient  to 

f 

match  the  transverse  observations. 

The  moment  tensors  in  the  near-field  from  this  inversion  are  given  in 
Figure  80  and  summarized  in  Table  2.  The  diagonal  elements  of  the  moment  ten¬ 
sor  are  no  longer  equal  in  magnitude  or  similar  in  time  function.  The  isotropic 
component  of  the  moment  tensor  is  a  relatively  simple  pulse  approximately  0.5 
to  0.67  seconds  in  duration.  The  magnitude  of  the  Isotropic  source  is  2.61 
times  larger  than  the  whole  record  inversion.  The  devlatoric  components  are 
more  complex  with  energy  spread  throughout  the  time  window.  The  ratio  of  the 
isotropic  to  the  maximum  devlatoric  components  is  only  1.14.  By  eliminating 
the  shear  waves,  the  fits  are  Improved,  the  isotropic  source  strength  increased, 
and  the  isotropic  to  devlatoric  ratio  decreased.  The  isotropic  source  spectrum 
is  given  in  Figure  81.  The  high  frequency  slope  is  3  with  a  corner  frequency 
at  2  Hz. 

A  summary  of  similar  P  wave  only  inversions  for  Handley  and  Pipkin  is 
given  in  Table  2.  As  found  for  Farm,  the  isotropic  source  Is  increased  with 
the  short  window  inversions,  while  the  Isotropic  to  devlatoric  ratio  decreases 
dramatically.  It  appears  that  the  shear  arrivals  are  very  important  In 
constraining  these  sources. 
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DISCUSSION/CONCLUSIONS 


Data  from  three  contained  nuclear  explosions  in  tuff  at  the  Nevada  Test 
Site,  Pahute  Mesa,  have  been  analyzed  using  the  moment  tensor  source 
characterization.  These  explosions  included  sources  just  above,  just  below, 
and  far  below  the  water  table.  An  average  velocity  model  has  been  developed 
in  order  to  complete  the  data  analysis.  This  model  utilizes  borehole  data, 
geological  information,  and  refraction  data.  The  velocity  structure  was 
further  refined  by  the  arrival  time  data  from  the  explosions.  The  most  impor¬ 
tant  features  of  the  model  are  the  comp>*ess1onal  velocity  jump  at  the  0.65  km 
water  table  followed  by  a  gradient  structure  (synthesized  by  2  layers)  down  to 
a  2.50  km  interface  that  acts  to  return  down  going  energy  to  the  free  surface. 

The  moment  tensor  source  formulation  has  been  used  to  fit  the  obser¬ 
vational  data  with  Green's  functions  from  the  just  described  velocity  model. 
Frequency  domain  inversions  for  the  complete  moment  tensor  and  just  its 
isotropic  component  were  done.  In  general,  the  vertical  and,  to  a  lesser 
degree,  the  radial  velocities  were  well  modeled  with  the  transverse  fits 
slightly  degraded.  The  moment  tensor  representation  can  be  used  to  adequately 
model  the  observational  data  in  the  near-field  (2-15  km)  from  underground 
explosions. 

Comparison  of  the  isotropic  moment  tensor  for  Pipkin,  Handley  and  Farm 
(Figure  82)  indicates  that  this  part  of  the  source  function  may  be  a  good 
measure  of  relative  explosive  yield.  The  measure  is  relative  If  the  propaga¬ 
tion  path  is  similar  for  the  explosions  under  study  and  it  can  become  absolute 
within  the  accuracy  of  Green's  functions  used  to  model  propagation  path 
effects . 

The  variety  of  Green’s  functions  used  in  this  study  show  shear  waves  which 
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average  5  to  IS  times  larger  than  the  explosion  P  waves.  In  order  to  obtain 

transverse  motions  comparable  to  radial  and  vertical  motions,  one  needs  a 

% 

t 

devlatorlc  source  which  Is  much  smaller  than  the  Isotropic.  In  the  moment 
tensor  Inversions,  these  phenomena  are  exactly  what  occurs.  In  fitting  the 
Farm,  Handley,  and  Pipkin  data  the  resulting  sources  are  dominated  by  the 
Isotropic  component  of  the  moment  tensor.  As  the  transverse  shear  waves  In 
the  Green's  functions  increased  relative  to  the  explosion  P  waves,  the  domi¬ 
nance  of  the  isotropic  moment  tensor  also  increased.  The  -olds  model  has  the 
largest  shear  waves  and  the  largest  ratio  of  isotropic  to  devlatorlc  source. 

The  variation  of  the  resulting  source  functions  between  the  three  explo¬ 
sions  and  with  changing  Green's  functions  was  summarized  In  Figures  76  and  77 
and  Table  2.  Of  all  the  source  characteristics,  the  corner  frequency,  high 
frequency  roll-off,  and  spectral  peak  are  least  variable.  The  source  phase 
Information  Is  somewhat  more  changeable  and  dependent  upon  the  Green's 
functions.  The  stability  of  the  spectral  information  In  these  few  Inversions 
suggests  exploration  of  spectral  characteristics  of  the  moment  tensor  for 
yield  and  discrimination  studies. 

The  Isotropic  time  functions  for  Farm,  Handley  and  Pipkin  all  suggest  the 
possibility  of  two  source  time  pulses  separated  by  1.0  to  1.5  seconds.  This 
observation  is  particularly  apparent  in  the  Farm  inversion.  The  physical 
interpretation  of  this  mathematical  representation  of  the  source  Is  open.  The 
secondary  source  could  be  due  to  nonlinear  processes  close-in  to  the 
explosion,  problems  In  our  propagation  model,  such  as  Improper  surface  wave 
calculation  or  data  preparation  problems.  These  areas  need  to  be  explored. 

The  Importance  of  shear  waves  In  controlling  the  isotropic  to  devlatorlc 
source  ratio  has  been  discussed  In  the  context  of  the  Green's  functions.  The 
inversions  utilizing  only  the  Initial  P  waves  from  the  observational  data 
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further  explored  this  question.  The  P  wave  only  inversions  yielded  greatly 
improved  fits,  but  the  isotropic  to  deviatoric  ratio  was  degraded.  It  is  felt 

f 

that  the  inclusion  of  shear  waves  in  the  studies  improves  our  constraints  on 
both  the  isotropic  and  deviatoric  components  of  the  source. 

A  preliminary  study  of  the  utility  of  the  moment  tensor  representation  in 
near-field  explosion  modeling  has  been  completed:  The  results  Indicate  the 
method  is  able  to  explain  the  observational  data- and  give  insight  into  rela¬ 
tive  and  possibly  absolute  yield  estimates.  The  puzzling  observation  of 
transverse  motions  as  large  as  radial  and  vertical  motions  may  be  partially 
explained  in  terms  of  the  strength  of  shear  waves  from  deviatoric  sources 
being  much  larger  than  the  isotropic  P  waves.  Such  models  allow  sources  which 
are  dominated  by  the  isotropic  component  and  have  transverse  motions  equal  to 
observed  radlals  and  verticals.  The  physical  Interpretation  of  the  source 
details  awaits  further  study.  This  Interpretation  should  take  advantage  of 
any  close-in  Instrumentation  in  the  nonlinear  ground  motion  regime.  Finally, 
the  stability  of  the  source  modulus  suggests  the  use  of  moment  tensor  spectral 
characteristics  in  discrimination  and  yield  studies. 
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ABSTRACT 


A  source  characterization  of  bermed  surface  bursts  utilizing  a  vertical 
and  horizontal  point  source  has  been  developed.  This  representation  Is  for¬ 
mulated  so  that  with  analytic  propagation  path  models  and  observational  data 
Inverse  techniques  can  be  applied  and  the  source  estimated.  The  inverse  pro¬ 
cedure  Is  applied  to  the  Pre-HYBRID  GUST  1-04,  05,  and  06  data  sets  (13.6  kg 
for  04,  13.6  kg  for  05,  and  39  kg  for  06  of  sand  bermed  C-4).  The  method  Is 
able  to  model  the  observational  data  explaining  the  R*l»®  decay  of  the  radial 
velocity  as  nearfield  body  waves  and  the  R~*5  decay  of  the  vertical  velocity 
as  Rayleigh  waves.  The  radial  source  function  Is  a  factor  of  2  larger  than 
the  vertical  for  all  tests.  Indicating  the  dominance  of  the  radial  force  In 
cratering.  Source  rise  time  correlates  well  with  the  time  the  sand  berm  con¬ 
tains  the  explosive  products  while  pulse  width  may  be  driven  by  nonlinear 
processes ,  such  as  spall.  The  sources  appear  to  cube  root  yield  scale. 
Stresses  predicted  using  the  determined  source  functions  Indicate  large  ver¬ 
tical  tensile  stresses  can  develop  from  the  Rayleigh  wave. 


INTRODUCTION 


In  an  Investigation  of  energy  coupling  from  surface  explosions,  the  Air 
Force  Weapons  Laboratory  has  detonated  a  number  of  single-  and  multiple-burst 
bermed  explosions.  The  surface  tangent  explosives  were  covered  with  berms  of 
sand  to  eliminate  the  alrblast  portion  of  loading,  leaving  only  the  energy 
directly  coupled  Into  the  ground  at  the  crater.  Alrblast  measurements  made 
during  these  tests  Indicate  the  bermlng  procedure  was  successful. 

It  Is  this  data  set  that  shall  be  explored  to  help  delineate  the  seismic 

source  characterl2at1on  that  Is  appropriate  for  the  cratering  process.  Much 

work  has  been  done  In  characterizing  contained  explosions,  attempting  to 

separate  source  and  propagation  path  effects  (Haskell,  1967;  McEvIlly  and 

Peppln,  1972;  Mueller  and  Murphy,  1971;  Werth  and  Herbst,  1963;  Hadley  and 

Helmberger,  1981;  Stump  and  Johnson,  1977,  1981).  Little  attention  has  been 
# 

paid  to  the  case  of  surface  bursts  (Ewing  and  Press,  1951;  Gupta  and 
Klssllnger,  1966).  For  ground  motion  prediction  purposes,  the  usual  approach 
to  the  surface  burst  problem  Is  to  take  the  data  of  contained  events,  correct 
for  the  geometric  effect  of  the  free  surface  on  the  propagation  path  and  then 
apply  the  results.  No  formal  consideration  Is  given  to  the  change  In  source 
characterization  as  the  explosion  approaches  the  free  surface.  To  first  order, 
the  contained  explosion  can  be  modeled  as  a  spherical  pressure  pulse  where  the 
bermed  surface  source  Is  modeled  as  a  point  force.  Elastic  responses  for 
these  two  sources  are  quite  different— the  contained  explosion  being  dominated 
by  P  waves,  while  the  point  force  at  the  free  surface  Is  rich  In  shear  and 
Rayleigh  waves  In  the  near  field. 

In  this  paper,  observational  data  from  these  bermed  surface  sources  are 
Investigated  taking  Into  account  the  difference  In  source  characterization 
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from  contained  and  bermed  surface  explosions.  Attempts  to  fit  the  obser- 
vatlonal  data  allows  the  refinement  of  the  surface  source  model. 

Finally,  this  study  should  shed  some  light  on  spall  mechanisms.  Spall,  as 
Indicated  by  -lg  dwells  on  vertical  accelerometers  followed  by  Impulsive 
rejoins  on  the  vertical  and  radial  accelerations,  were  observed  on  all  the 
bermed  surface  explosions  (Stump  and  Relnke,  1982).  Since  the  source  Is  at 
the  free  surface,  the  spall  or  tensile  failure  of  the  material  cannot  be 
explained  In  terms  of  a  compresslonal  wave  leaving  the  source,  traveling  to 
the  free  surface,  and  reflecting  as  a  tensile  wave  (Elsler  and  Chilton,  1964). 
Both  the  source  characterization  and  the  propagation  path  effects  must  be  well 
accounted  for  In  order  to  adequately  propose  a  mechanism  for  tensile  failure. 
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SOURCE  REPRESENTATION 

Assuming  a  source  can  be  represented  as  either  a  body  force.  Initial 
condition,  or  boundary  condition  then  a  general  representation  can  be  written 


UfcU'.t')  Is  the  ground  motion  In  the  k  direction  at  (x',t'),  V0  Is  the  source 
volume,  Gklfx'.t'jx.t)  is  the  displacement  in  the  k  direction  at  (x',t')  due  to 
a  delta  function  force  In  the  1  direction  at  (x,t),  and  ff  Is  the  equivalent 
body  force  in  the  1  direction  at  (x,t.).  The  summation  convention  is  assumed. 

In  the  case  of  a  contained  explosion,  assuming  the  wavelengths  of  Interest 
are  long  compared  to  the  characteristic  source  length,  equation  (1)  becomes: 


Assuming  a  bermed  surface  source  can  be  modeled  as  a  vertical  and  horizon¬ 


tal  point  force,  then  the  representation  becomes: 

CO 

U*A*\*)  *  (3) 

The  1  direction  has  been  chosen  radial  away  from  the  source  and  3  direction  Is 


positive  down. 

Transforming  equations  (2)  and  (3)  Into  the  frequency  domain  eliminates 
the  time  Integral  and  allows  one  to  write: 

U  ■  GF  (4) 


U  Is  the  observed  ground  motion,  G  Is  the  propagation  path  effect ,  and  F  Is 
the  source  characterization.  For  the  surface  source,  F  consists  of  two  terms. 


one  for  the  vertical  and  one  for  the  horizontal  force.  In  the  case  of  the 
buried  explosion,  F  consists  of  six  unique  terms. 


Equation  4  is  written  for  each  frequency  so  that  if  there  are  256  initial 
time  points,  then  the  equation  Is  written  128  times  (real  and  Imaginary)  for 
each  frequency  point.  For  each  frequency,  the  U  matrix  Is  made  up  of  the  real 
and  imaginary  components  of  the  Fourier  transformed  ground  motion  data  and  Is 
2n  x  1  in  dimension  (where  n  Is  the  number  of  Individual  observations);  the  G 
matrix  Is  composed  of  the  real  and  Imaginary  parts  of  the  Fourier  transformed 
Green's  functions  and  is  2n  x  4  In  dimension;  and  the  F  matrix  Is  4  x  1  and 
consists  of  the  real  and  Imaginary  parts  of  Fj  and  F3  (equation  3). 

In  the  study,  U  are  the  observed  velocities  In  the  frequency  domain  from 
the  PHG  tests  and  G  the  propagation  path  effects  determined  analytically  for 
an  elastic  half-space.  The  only  unknowns  are  F,  the  source  characterization. 
Calculating  the  inverse  of  the  matrix  G  and  multiplying  It  by  U  allows  one  to 
determine  F. 

F  -  G^U  (5) 

The  problem  Is  underdeterml ned  when  there  are  less  than  two  seismograms,  exact 
when  there  are  two  seismograms,  and  may  be  over determined  when  there  are  more 
than  two  observations.  The  Inversions  In  this  study  will  use  a  maximum  of  ten 
seismograms  and  a  minimum  of  six. 

In  solving  Equation  4  and  writing  Equation  5,  singular  value  decomposition 
of  6  Is  used  and  a  generalized  Inverse  Is  calculated  (Lanczos,  1961).  G  Is 
decomposed  as: 

G  -  wqvT  (6) 

Where  W  Is  a  2nxp  semi -orthogonal  matrix  consisting  of  the  eigenvectors  asso¬ 
ciated  with  the  nonzero  eigenvalues  of  GGT,  V  Is  the  2  x  p  semi orthogonal 
matrix  consisting  of  similar  eigenvectors  of  GTG,  and  Q  Is  a  pxp  matrix  whose 
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diagonal  members  are  the  positive  square  roots  of  the  largest  p  nonzero  eigen- 
values  of  G^G. 

The  generalized  inverse  of  G  becomes: 

8-1  .  VQ-1  WT  (7) 

This  formulation  always  exists,  since  Q-*  can  always  be  ca^rlated  with 
(eigenvalues  of  G^G)'^  down  the  matrix  diagonal-  For  the  case  when  p  is  less 
than  the  dimmenslon  of  G^G,  the  Inverse  ignores  the  eigenvectors  associated 
with  the  zero  eigenvalues  and  thus  may  not  give  a  unique  solution  to  the 
problem.  The  inverse  Is  unable  to  resolve  the  solution  In  the  directions  of 
the  eigenvectors  associated  with  the  zero  eigenvalues.  The  resolution  matrix 
can  be  calculated  as  G"lG  «  WT.  The  eigenvalues  and  eigenvectors  of  this 
matrix  tell  one  If  the  problem  Is  well  enough  posed  to  allow  the  unique  recov¬ 
ery  of  the  source  characterization. 

Once  F  has  been  determined  using  Equation  5,  predicted  ground  motions 
(acceleration,  velocity,  displacement)  can  be  estimated  by  substituting  F  Into 
Equation  3  with  the  appropriate  Green's  functions.  The  stress  tensor  can  4  ;en 
be  calculated  as  follows: 


j ‘S'.-t')  *  P-  feu  +  Gj,,L  I*',* 

’  r  (8) 

n-  (&■  .  u'.v +• 


It  should  be  emphasized  that  once  F  has  been  determined  from  Equation  5  the 
stresses  and  motions  can  be  predicted  anywhere  In  the  material  using  Equations 
8  and  3. 
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GREEN'S  FUNCTIONS 


The  propagation  path  effects  In  this  study  are  modeled,  as  an  elastic 
half-space.  The  motivation  for  this  choice  Is  threefold:  (1)  The  radial  and 
vertical  velocity  waveforms  from  all  of  the  explosions  show  relatively  simple 
waveforms  with  the  radial  component  of  motion  dominated  by  body  waves  and  the 
vertical  component  primarily  surface  waves  (Figure  1),  (2)  the  body  wave 
decay  rate  of  R”**®  Indicates  the  strong  influence  of  near-field  components  of 
motion,  and  (3)  the  half-space  Green's  functions  allow  the  analytic  calculation 
of  the  resulting  stress  tensor.  The  half-space  model  Is  further  supported  by 
geophysical  work  (Stump  and  Relnke,  1982)  and  the  first  arrival  time  data  at 
the  test  site.  There  Is  a  hint  in  the  data  of  the  presence  of  a  layer  over  a 
half-space  structure.  The  first  arrival  time  data  Indicate  a  possible  break 
over  point  at  about  10  m  while  the  first  motion  of  the  vertical  motion  Is  up. 
This  motion  is  small  compared  to  the  predominantly  radial  body  waves  at  these 
distances.  Our  philosophy  Is  to  accept  the  simplest  propagation  and  source 
model  that  explains  the  observational  data  base. 

The  wave  front  diagram  for  a  surface  source  and  buried  receiver  In  an 
elastic  half-spac.e  Is  given  In  Figure  3.  The  first  arrival  Is  primarily  a 
radial  P  wave.  The  second  arrival  Is  diffracted  PS.  This  wave  leaves  the 
source  as  a  P  wave  traveling  along  the  free  surface  and  then  refracts  Into  the 
material  as  an  S  wave  to  the  receiver.  Following  these  arrivals  are  the  shear 
and  Rayleigh  waves. 

The  Green's  functions  for  a  surface  point  force  source  and  a  burled 
receiver  are  given  In  Figure  4.  In  this  model,  the  P  wave  velocity  Is  366 
m/sec  with  a  Poisson's  ratio  of  0.10.  The  seismic  velocities  were  determined 
from  refraction  studies  In  the  area  (Stump  and  Relnke  1982).  The  delta  func¬ 
tion  displacement  responses  for  a  radial  (fj),  transverse  (fg),  and  vertical 
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(f3)  point  force  source  are  Illustrated.  The  P,  diffracted  PS,  S,  and 
Rayleigh  arrivals  are  noted.  It  Is  the  Green's  functions  due  to  fj  and 
f3  which  will  be  substituted  Into  G  of  Equation  4  for  determination  of  F. 

The  effect  of  reducing  the  shear  wave  velocity  while  holding  the  P  wave 
velocity  constant  is  given  In  Figure  5.  These  Green's  functions  are  the 
radial  and  vertical  displacements  due  to  a  radial  force  localized  in  space  and 
time.  As  the  shear  wave  velocity  increases,  the  amplitude  of  the  shear  and 
Rayleigh  waves  increases.  At  large  Poisson's  ratio,  an  additional  arrival, 
the  so-called  Leaking  mode  P,  is  seen  on  the  radial  motion  between  the  P  and  S 
wave.  This  phenomena  was  discussed  by  Gilbert  and  Knopoff  (1961),  Gilbert,  et 
al.  (1962),  and  Chapman  (1972).  For  a  complete  discussion  of  this  arrival  In 
the  context  of  both  motion  and  stress  calculations,  see  Stump  (1982). 
Inversions  utilizing  Green's  functions  with  Poisson's  ratios  0.10  to  0.45  will 
be  accomplished  to  Investigate  the  trade-off  between  the  source  and  propaga¬ 
tion  effects. 
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DATA  ANALYSIS 


The  focus  of  the  source  characterization  In  this  study  was  a  series  of 

f 

bertned  explosions.  The  test  series,  known  as  Pre-HYBRID  GUST  (PHG),  contained 
three  bermed  sources.  These  sources,  PHGI-04,  PHGI-05,  and  PHGI-06,  consisted 
of  13.6  kg,  13.6  kg,  and  39  kg  of  C-4  explosive  under  a  sand  berm.  Figure  6 
summarizes  the  explosive  sources  and  sand  berm  size  which  changed  with 
experiment. 

The  three  sources  were  Instrumented  with  radial  and  vertical  accelerome¬ 
ters  in  the  3  to  36  m  range.  The  data  was  recorded  analog  on  high  speed  FM 
tapes.  These  records  were  then  digitized  at  2,000  samples  per  second  with  an 
antialias  filter  at  400  hz  (5  Pole  Butterworth  filter).  Figure  7  Illustrates 
a  sample  spectra  (solid  line)  and  noise  estimate  (dotted  line)  for  the  radial 
gage  at  6.55  m  from  PHGI-06.  Good  signal  to  noise  ratios  were  found  out  to 
400  hz.  Both  the  acceleration  and  integrated  velocity  records  are  given.  In 
the  inversions,  the  velocity  data  Is  used  with  a  correction  for  a  ramp.  This 
record  Is  characteristic  of  data  quality  for  all  experiments. 

1.  PHGI-06 

The  first  data  set  to  which  the  Inversion  scheme  was  applied  was  PHGI-06. 
This  experiment  was  the  biggest  of  the  three  studied  and  consisted  of  39  kg  of 
the  explosive  C-4.  Vertical  and  radial  data  at  the  6.55,  8.50,  10.85,  14.60, 
and  18.35-meter  ranges  (1.50  m  depth)  were  used.  The  Instrumental  array  and 
data  characteristics  are  summarized  In  Figure  8.  Prior  to  the  Inversion  and 
after  antialias  filtering,  the  velocity  data  was  decimated  to  1,000  samples 
per  second. 

All  ten  components  of  data  were  used  In  the  Inversion  to  determine  the  two 
equivalent  body  forces  fj  and  f3>  Both  the  data  and  the  half-space  Green's 
functions  were  Fourier  transformed  prior  to  substitution  into  Equation  4.  The 
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U  matrix  for  each  frequency  became  20x1,  while  the  G  matrix  was  20x4.  Since 
256  msec  of  data  was  analyzed,  the  20x4  G  matrix  was  Inverted  128  times. 

After  f}  and  f3  were  determined  using  Equation  5,  F  and  G  were  substituted 
back  Into  Equation  4  to  yield  a  predicted  U  for  comparison  to  the  obser¬ 
vational  data.  The  similarity  of  the  observed  U  to  the  calculated  U  Is  a 
measure  of  the  adequacy  of  the  source  model  and  Inversion.  A  summary  of  peak 
velocity  amplitudes  in  the  observed  and  calculated  waveforms  Is  given  In 
Figure  9.  The  radial  data  and  fits  decay  as  R“l*6,  while  the  verticals  decay 
as  The  vertical  data  are  particularly  well  matched,  while  errors  as 

large  as  a  factor  of  two  can  be  seen  in  the  radials.  Taking  the  data  set  as  a 
whole,  both  the  radial  and  vertical  peaks  are  well  fit. 

Now  focusing  on  the  entire  wave  train,  the  observational  and  predicted 
waveforms  from  the  Inversions  are  illustrated  in  Figures  10,  11,  and  12.  In 
each  plot,  the  observed  and  calculated  waveforms  are  plotted  to  the  same 
scales.  The  absolute  scale  for  each  plot  may  change  from  gage  to  gage  with 
the  peak  velocity  for  each  gage  Indicated  In  plots.  A  total  of  256  msec  of 
data  were  fit  and  plotted.  The  sign  convention  for  the  data  Is  positive 
radial  away  and  vertical  down.  Tne  fits  of  the  observational  waveforms  are 
well  explained  by  the  source  model.  The  mean  correlation  coefficient  between 
the  observed  and  calculated  velocities  Is  0.84.  The  radial  components  of 
velocity  are  dominated  by  the  near  field  body  wave  yielding  the  R-l*6  spatial 
decay.  The  vertical  component  of  velocity  is  explained  by  the  Rayleigh  wave 
arrival  giving  R"0*5  decay.  The  most  glaring  deficiency  In  the  fits  Is  the 
inability  to  model  the  small  Initial  upward  pulse  on  the  vertical  component  of 
motion. 
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As  Figure  9  Illustrates,  the  radial  component  of  velocity  Is  largest  at 
distances  less  than  7  m  with  the  vertical  component  large?  at  greater 
distances.  There  Is  some  hint  at  a  change  In  decay  rate  for  the  radial  motion 
at  approximately  18  m.  At  this  distance,  the  Rayleigh  wave  on  the  radial  com¬ 
ponent  Is  approaching  the  same  size  as  the  body  wave.  One  can  hypothesize 
that  at  ranges  greater  than  18  m  the  Rayleigh  wave  will  dominate  the  radial 
motion,  and  spatial  decay  will  be  close  to  R”»S. 

The  source  time  functions  and  their  Integrals  for  the  radial  (f^)  and  ver¬ 
tical  (f3)  point  forces  determined  from  the  PHGI-06  inversions  are  given  in 
Figure  13.  As  indicated  by  the  relative  source  strengths,  the  radial  force  is 
a  factor  of  two  greater  than  the  vertical.  The  net  impulse  of  the  sources  is 
down  and  away  from  ground  zero.  There  is  some  problem  with  causality  In  the 
smaller  vertical  time  function. 

In  the  time  domain,  the  source  rise  time  is  18  msec  with  a  pulse  width  of 
70  msec.  Ake  (1980)  in  an  analysis  of  the  PHGI-06  films  determined  the  time 
to  venting  of  the  berm  to  be  17  msec.  Figure  15  illustrates  the  berm  growth 
as  a  function  of  time  with  the  explosive  by  products  breaking  out  of  the  berm 
at  17  msec.  It  appears  that  the  source  rise  time  Is  controlled  by  the  time 
that  the  sand  berm  contains  the  explosion. 

The  corresponding  source  spectra  are  given  in  Figure  14.  The  factor  of 
two  difference  in  source  strength  between  fj  and  f3  can  be  seen  in  the 
spectra.  The  corner  between  10  and  20  hz  Is  representative  of  the  pulse  width 
observed  in  the  data.  The  high  frequency  roll-off  of  the  source  Indicates  a 
slope  of  between  2  and  3. 

The  least  known  parameter  In  the  Inversion  scheme  Is  the  shear  wave  velo¬ 
city  of  the  test  bed.  In  the  Initial  inversions,  a  Poisson's  ratio  of  0.10 
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was  assumed.  To  Investigate  the  stability  of  the  Inversion  procedure  and  the 
effect  of  varying  the  shear  wave  velocity,  two  additional  runs  were  made  with 
Poisson's  ratios  of  0.25  and  0.40.  The  effect  of  varying  Poisson's  ratio  on 
the  Green's  functions  was  shown  In  Figure  5.  Not  only  does  the  P  to  Rayleigh 
wave  amplitude  decrease  as  Poisson's  ratio  Increases,  but  the  characteristic 
frequency  of  the  whole  waveform  also  decreases.  The  relative  sizes  of  fj  and 
f3  for  the  Inversions  using  the  three  sets  of  Green's  functions  are  given  In 
Figure  16.  The  two  to  one  ratio  of  fj  to  f3  remains  the  same  for  all 
inversions.  The  relative  source  strength  changed  by  a  factor  of  three  when 
Poisson's  ratio  changed  from  0.10  to  0.25.  This  change  was  only  a  20  percent 
variation  in  shear  wave  velocity.  Source  coupling  for  this  characterization 
is  a  strong  function  of  Poisson's  ratio.  The  time  functions  of  fi  and  for 
the  three  Poisson's  ratios  are  given  In  Figures  17  and  18.  Since  the  Green's 
functions  decrease  In  frequency  as  Poisson's  ratio  Increases,  the  resulting 
source  time  functions  Increase  In  frequency.  The  effect  of  the  changing  Green’s 
functions  on  the  fits  to  the  observational  data  Is  summarized  for  the  vertical 
velocity  at  10.85  m  In  Figure  19.  The  fits  for  Poisson's  ratio  of  0.10  and 
0.25  are  Indistinguishable,  while  the  fit  for  0.40  Is  greatly  degraded.  Using 
the  fits  as  a  guide  to  the  adequacy  of  the  model,  one  cannot  distinguish  be¬ 
tween  a  Poisson's  ratio  of  0.10  and  0.25.  This  conclusion  leaves  a  relative 
error  in  source  coupling  of  3.  It  Is  felt  that  the  20  percent  change  In  shear 
wave  velocity  encompassed  by  these  models  Is  difficult  to  resolve  through 
geophysical  exploration  techniques. 

2.  PHGI-05 

The  second  data  set  to  which  the  source  characterization  scheme  was 
applied  was  PHGI-05  (Figure  6).  This  experiment  was  identical  to  PHGI-06, 
except  the  yield  was  reduced  from  39  kg  to  13.6  kg  of  C-4. 
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The  Instrumentation  was  somewhat  more  scarce  In  this  test,  so  only  three  gage 
locations  were  used  In  the  Inversion— 6.55,  8.50,  and  18.35  meters  (Figure 
20).  The  data  analysis  was  conducted  with  Green's  functions  for  a  Poisson's 
ratio  of  0.10.  This  final  propagation  path  model  was  chosen  because  of  Its 
agreement  with  geophysical  exploration  work  at  the  test  site  and  the  success 
of  fits  to  the  PHGI-06  data.  As  the  Poisson's  ratio  tests  on  the  PHGI-06  data 
show,  Poisson's  ratio  can  go  as  high  as  0.25. 

Using  the  same  format  as  In  the  PHGI-06  analysis,  the  observed  and  calcu¬ 
lated  seismograms  are  given  In  Figures  21  and  22.  As  In  the  previous  results, 
the  fits  are  quite  good. 

The  source  functions  and  their  spectra  are  In  Figures  23  and  24  respec¬ 
tively.  The  rise  time  of  the  source  Is  again  18  msec  which  compares  to  Ake's 
(1980)  observation  of  berm  venting  at  19.5  msec  (Figure  25).  The  2.5  msec 
difference  In  vent  time  between  PHGI-05  and  PHGI-06  Is  too  small  to  resolve  by 
this  procedure.  The  pulse  width  of  the  source  Is  approximately  64  msec, 
slightly  smaller  than  the  70  msec  observed  In  PHGI-06. 

Since  the  berm  size  for  PHGI-05  and  PHGI-06  Is  similar,  one  can  begin  to 
study  the  effect  of  source  yield  on  coupling.  Comparison  of  Inversions  uti¬ 
lizing  Poisson's  ratio  of  0.10  Is  made  In  Figure  26.  The  dotted  line  In  the 
figure  is  that  predicted  by  cube  root  yield  scaling.  The  limited  source 
strength  data  fall  close  to  this  line. 

3.  PHGI-04 

The  final  data  set  analyzed  In  this  study  was  PHGI-04  (Figure  6).  The 
experiment  was  Identical  to  PHGI-05,  but  the  berm  was  reduced  In  height  and 
width  by  20  percent.  Again,  three  gage  locations  at  6.55,  8.5  and  18.35 
meters  were  used  In  the  Inversion  (radial  and  vertical  velocity  at  each 
location).  As  In  the  PHGI-05  Inversion,  the  Green's  functions  were  for  a 


Poisson’s  ratio  of  0.10. 

The  observed  and  calculated  seismograms  are  given  In  Figures  27  and  28. 

As  In  the  other  tests,  the  fits  are  good.  The  source  time  functions  and  their 
spectra  are  given  In  Figures  29  and  30  respectively.  From  the  films.  Figure 
31,  the  time  to  venting  was  12.5  msec  (Ake,  1980).  This  time  compares  to 
11-12  msec  rise  time  on  the  source  function.  As  In  the  previous  inversions, 
the  fj  to  f3  ratio  remained  approximately  2. 


STRESS  PREDICTIONS 


Using  the  sources  determined  from  the  observations,  one  can  predict  the 
stress  tensor  by  plugging  fi(t)  and  f3(t)  Into  Equation  8.'  One  can  then  pre¬ 
dict  the  stresses  any  where  within  the  test  bed  (Stump,  1981). 

Using  the  source  from  PHGI-06  Inversions,  the  stress  tensor  at  a  distance 
of  10.85  m  and  a  depth  of  1.37  m  were  predicted  (Figure  32).  In  these  plots, 
down  Is  compression  and  up  Is  tension.  The^u  stress  Is  the  largest  with  a 
compressive  pulse  controlled  by  the  body  wave  arrival  (the  1  direction  Is 
radial  away  from  the  source).  A  tensile  hoop  stress  can  be  seen  on  the 622 
component.  The 633  component  yields  a  large  tensile  stress  In  the  vertical 
direction  as  a  result  of  the  Rayleigh  wave.  This  tensile  stress  may  be  a 
mechanism  leading  to  spall  from  surface  explosions.  The  stress  tensor  Is 
decomposed  Into  Its  devlatorlc  and  Isotropic  components  In  Figure  33. 

Similar  stress  calculations  were  done  at  a  variety  of  distances  and  depths 
from  the  source.  A  summary  of  the  maximum  vertical  tensile  stress  from  these 
calculations  Is  given  In  Figures  19  and  20.  The  decay  of  the  stress  with 
range  from  the  source  for  a  depth  of  1.50  m  Is  given  In  Figure  34.  In  Figure 
35,  the  stresses  are  found  to  Increase  with  depth  for  the  10.85  m  range.  This 
apparent  Increase  of  stress  with  depth  Is  a  result  of  the  zero  stress  boundary 
condition  at  the  free  surface.  Since  the  tensile  stress  Is  a  result  of 
Rayleigh  waves,  the  33  term  will  decay  with  depth.  None  of  the  calculations 
shown  were  deep  enough  to  eliminate  the  free  surface  effect  (at  4  m,  one  Is 
only  1/5  of  a  wavelength  from  the  free  surface). 
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CONCLUSIONS 


The  waveforms  from  the  PHG  single-burst  bermed  explosions  can  be 
Interpreted  as  being  dominated  by  nearfleld  body  wav&s,  R*^*6,  on  the  radial 
components  and  surface  waves,  R"*5,  on  the  vertical  velocity  components. 

Using  this  Interpretation  and  a  vertical  and  horizontal  point  source  model, 
one  Is  able  to  fit  the  data  with  an  average  correlation  coefficient  on  PHGI-06 
of  0.84. 

These  source  Inversions  Indicate  the  source  coupling  Is  sensitive  to 
Poisson's  ratio.  A  20  percent  change  In  shear  wave  velocity  (Poisson's 
ratio  0.10  to  0.25)  resulted  In  a  factor  of  3  change  In  coupling.  This 
variation  In  shear  properties  Is  what  we  expect  In  the  variability  of  our 
shear  properties  from  geophysical  exploration.  We,  therefore,  conclude  that 
the  factor  of  3  error  In  coupling  estimate  Is  a  lower  bound  for  these  sources. 

A  visual  comparison  of  the  source  time  functions  from  the  three  bermed 
experiments  Indicates  great  similarities  (Figures  36  and  37).  PHGI-06  and 
PHGI-05  have  nearly  Identical  time  functions,  although  their  relative 
strengths  Indicate  the  decrease  In  absolute  yield  from  39  kg  to  13.6  kg 
(PHGI-06  to  PHGI-05).  The  smaller  bermed  explosion,  PHGI-04,  shows  a  higher 
frequency  source  than  the  two  previous.  In  all  cases,  the  rise  time  of  the 
source  agrees  within  one  millisecond  of  the  time  to  break  out  of  explosive 
products  from  the  sand  berm.  The  pulse  widths  of  the  sources  are  much  longer 
In  duration  than  the  rise  time  and  may  correlate  with  other  nonlinear 
processes,  such  as  spall.  The  ratio  of  radial  to  vertical  forces  of  2  to  1 
remained  constant  for  all  sources  using  a  variety  of  Green's  functions.  It  Is 
felt  that  this  ratio  represents  the  dominance  of  the  radial  force  In 
cratering.  The  mathematical  representation  of  these  sources  has  been  well 
correlated  with  physical  phenomenology  near  the  explosion. 
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Using  the  source  function  developed  from  the  observational  ground  motion 
data,  the  predicted  stress  tensor  has  been  given.  The  stresses  yield  tensile 
components  In  the  vertical  direction  resulting  from  the  Rayleigh  wave  arrival. 
This  tensile  stress  Increases  with  depth  over  the  0-4  m  depth  range  due  to  the 
free  surface  effect.  A  mechanism  for  tensile  failure  of  a  material  due  to 
surface  bursts  has  been  Identified. 


16 


BIBLIOGRAPHY 


Ake,  J.P.,  1980.  Ground  Motion  In  a  Non-AIrblast  Environment  Letter  Report, 
Engineering  Research  Institute,  University  of  New  Mexico. 

Babcock,  S.,  1980.  Pre-HYBRID  GUST  Phase  I  Quick  Look  Report,  Engineering 
Research  Institute,  University  of  New  Mexico,  CERF-AG-26. 

Gupta,  I.N.  and  C.  Klssllnger,  1966.  Radiation  of  Body  Waves  from 
Near-Surface  Explosive  Sources,  Geophysics  31,  1057 

Hadley,  O.M.  and  O.V.  Helmberger,  1981.  Seismic  Source  Functions  and 

Attenuation  from  Local  and  Teleselsmlc  Observations  of  the  NTS  Event  Jorum 
and  Handley,  Bull.  Seism.  Soc.  Am.  71. 

Haskell,  N.A.,  1967.  Analytic  Approximation  for  the  Elastic  Radiation  From  a 
Contained  Underground  Explosion,  J.  Geophys.  Res.  72.  2583-2587. 

McEvllly,  T.V.  and  W.A.  Peppln,  1972.  Source  Characteristics  of  Earthquakes, 
Explosions,  and  After  Events,  Geophys.  J.  Roy.  Astr.  Soc.  31.  67-82. 

Mueller,  R.A.  and  J.R.  Murphy,  1971.  Seismic  Characteristics  of  Underground 
Nuclear  Detonations:  Part  I,  Seismic  Scaling  Law  of  Underground  Nuclear 
Detonations,  Bull.  Seism.  Soc.  Am.  61,  1675. 

Press,  F.  and  M.  Ewing,  1951.  Ground  Roll  Coupling  to  Atmospheric 
Congressional  Waves,  Geophysics  16.  416-430. 

Stump,  B.W.  and  L.R.  Johnson,  1977.  The  Determination  of  Source  Properties  by 
the  Linear  Inversion  of  Seismograms,  Bull.  Seism.  Soc.  Am.  67,  1489-1502. 

Stunp,  B.W.  and  L.R.  Johnson,  1981.  The  Effect  of  Green's  Functions  of  the 
Determination  of  Source  Mechanisms  by  the  Linear  Inversion  of  Seismograms, 
Identification  of  Seismic  Sources  -  Earthquake  or  Underground  Explosion. 
Edited  by  Dr  Eysteln  Husebye  and  S.  Mykkehvelt,  D.  Reldel  Publishing  Co., 
Dordrecht-Hol 1  and. 

Stump,  B.W.,  1982.  Stress  Waves  In  a  Three-Dimensional  Elastic  Half-Space: 
Single  and  Multiple  Surface  Sources,  In  preparation. 

Werth  and  Herbst,  1963.  Comparison  of  Amplitudes  of  Seismic  Waves  from 
Nuclear  Explosions  in  Four  Mediums,  J.  Geophys.  Res.  68,  1463-1475. 


17 


I  •  $lnV/**> 


Figure  3 


DISPLACEMENTS -6(1) 


0 


100 


200  h  •  1.37 
msec 


Figure  5 


20 


SPECTRUM 


22 


FREQUENCY  ( HZ  I  Figure  7 


PHG 


3  vertical  (down) 


PH6  1-06  FITS 
polsson's  ratio  0.10 
R  ■  6.55m 


.I64m/s 
observed  J 


radial 


jc&lculated 


T 


I 


calc 


j 


too 


200  msec 


Fi  gure  10 


PHG  I  -  06 
Berm  Growth 


msec 

msec 

msec 

msec 


5.64m 


Figure  15 


RELATIVE  SOURCE 
STRENGTH 


EFFECT  OF  POISSON’S  RATIO  ON 
SOURCE  STRENGTH  , 


4 

.1 


.2  3 

POISSON'S  RATIO 


Figure  16 
31 


RATIO  ON  SOURCE 


EFFECT  OF  POISSON'S  RATIO  ON  FITS 
PHG  1-06,  R- 10.85m  H-I5m  V 


.067m/s 


.064  m/s 


,064m/s 


.037m/s 


i _ i _ i _ j _ i _ i 

0  100  200msec 


Figure  19 


34 


PHG  1-05 


35 


3  vertical  (down) 


PHG  I  -  OS  PITS 
polsson's  ratio  O.IO 
R*  6.39m 


.096  m/s  T 

J&s  1 


radial 


t 

away 


CALC 


1 


Royleigb 


.102  m/s 
OBS_ 


I 


1 

down 


vortical 


PHG  1-05  SOURCE  FUNCTIONS 

POISSON'S  RATIO- 0.1 


—  I  ■  —  -  1  1  *  ** 

100  200  msec 

Figure  23 


38 


SOURCE  SPECTRA 


ncauoer'tfJ 


PHG  I  -  05 
Berm  Growth 


RELATIVE  SOURCE 
STRENGTH 


SOURCE  SCALING 


10  20  30  40 

YIELO  (kg) 


Figure  26 


41 


PHGI-04  FITS 


poisson *s  rotto  O.IO 
R  -  6.55m 

— 


0.07 2 m/s 
—  obs 


0.046m/s 
—  calc 


0.093  m/s 
—  obs 


0.057  m/s 
—  calc 


i _ ■  i  i  » _ i 

0  100  200  msec 


Figure  27 
42 


UIQ£'8I 


< 

cc 

o** 
UJO 
CL  | 
(fl¬ 
ed 
UIX 
O  0- 


SnTiQOw 


PHG  I  -  04 
Berm  Growth 


DEVIATORIC  STRESS 


COMPRESSION 


PHG  1-06  PRED. 


RANGE(m) 

figure  34 


49 


PH 6  1-06  STRESS  PRED, 


.5  + 


54+ 

</> 

jo  < 
m  m 

O 

of** 

W 

s*r 

$2 


m 


.1+ 


< - 1 - 1 - h 

12  3  4 

DEPTH  (m) 


Figure  35 
50 


SOURCE  SUMMARY  PHG  1-04,05,06 


Figure  37 


CONTAINED  HIGH-EXPLOSIVE 
ALLUVIUM  TEST 
DATA  REPORT 


BRIAN  W.  STUMP 
ROBERT  E.  REINKE 
AFWL/NTES 

KIRTLAND  AFB,  NEW  MEXICO  87117 


INTRODUCTION 


During  August  of  1981,  the  Air  Force  Weapons  Laboratory,  In. con junction 

f 

with  personnel  from  the  University  of  California  at  Berkeley,  Instrumented  and 
detonated  a  chemical  explosion  of  253  pounds  burled  at  a  contained  depth  of 
11.5  meters.  The  explosive  used  was  a  TNT  sphere  20%  inches  In  diameter.  The 
test  was  fired  in  alluvium. 

The  purpose  of  this  test  was  fourfold:  (1)  We  wished  to  conduct  a  highly 
instrumented  contained  explosive  shot  In  a  relatively  simple,  known  geology. 
Such  an  experiment  would  allow  a  careful  characterization  of  the  explosive 
source  from  a  seismic  point  of  view.  (2)  The  fully  contained  test,  in  con¬ 
junction  with  a  variety  of  depth  of  burst  tests  (already  completed)  at  the 
same  test  site,  would  allow  one  to  assess  the  importance  of  burial  depth  on 
source  characterization  and  relative  coupling.  (3)  The  fully  contained  che¬ 
mical  explosion  (HE),  in  conjunction  with  similar  nuclear  explosion  (NE)  data 
sets  from  the  Nevada  Test  Site,  would  allow  one  to  study  the  differences  of  NE 
and  HE  sources.  In  particular,  we  hoped  to  study  the  applicability  of  various 
scaling  relations  trying  to  develop  separate  procedures  for  source  and  propa¬ 
gation  path.  (4)  Spall  or  tensile  failure  of  the  soil  as  identified  by  -lg 
dwells  on  vertical  accelerometers  followed  by  impulsive  rejoins  on  the  ver¬ 
tical  and  radial  accelerations  was  identified  on  records  from  a  number  of  sur¬ 
face  and  near  surface  bursts.  The  mechanism  of  spall  for  these  tests  could 
not  be  explained  In  terms  of  a  compressive  wave  leaving  the  source  reflecting 
as  a  tensile  wave  at  the  free  surface  and  then  failing  the  material.  In  order 
to  tie  the  surface  explosion  spall  data  set  to  the  contained  nuclear  data,  a 
contained  chemical  explosion  was  suggested. 

We  intend  to  include  no  analysis  in  this  report.  The  purpose  of  this  work 
is  to  outline  the  experiment  and  present  the  observed  data. 
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TEST  OUTLINE 


In  order  to  begin  the  resolution  of  the  four  questions  already  outlined,  a 
contained  2S3-pound  (TNT)  chemical  explosive  test  was  designed.  The  test  site 
was  chosen  to  supplement  a  series  of  depth  of  burst  studies  already  completed. 
The  site  was  plane  layered  alluvium,  close  to  the  Air  Force  Weapons  Laboratory 
In  Albuquerque,  New  Mexico. 

The  Instrumentation  array  for  the  test  consisted  of  three  distinct  sets. 
The  first  part  of  the  instrumentation  was  made  up  of  close-in  accelerometers 
mainly  distributed  along  one  radial  (Figure  1).  The  primary  purpose  of  these 
gages  was  to  yield  good  spatial  control  on  the  nonlinear  processes  close  to 
the  source  and  then  transition  Into  the  near-field  part  of  the  Instrumentation 
which  will  be  used  for  the  seismic  characterization  of  the  source.  We  were 
particularly  interested  in  the  timing  and  propagation  direction  of  the  spall 
phase,  so  six  2-component  gages  were  placed  at  a  radius  of  1  m  from  GZ  in  a 
vertical  array.  In  order  to  maintain  some  azimuthal  control  close-in  to  the 
source,  two  additional  azimuths  were  instrumented  at  the  3-meter  range  (Figure 
2).  All  these  gages,  which  were  later  digitized,  were  recorded  on  analog 
tape.  The  FM  tape  drives  are  run  at  a  high  rate  of  speed,  since  the  shot 
time  is  known  and  the  signals  quite  short.  As  a  result,  good  signal  to  noise 
ratios  are  obtainable,  assuming  good  preshot  ground  motion  estimates  are  made. 
The  data  from  this  test  were  digitized  at  5000  samples  per  second  with  a 
5-pole  butterworth  filter  at  1250  Hz.  A  total  of  48  gages  were  fielded. 

The  so-called  near-field  array  was  fielded  with  the  help  of  Lane  Johnson 
and  Tom  McEvIlly  from  the  University  of  California  at  Berkeley.  The  primary 
design  of  this  part  of  the  Instrumentation  was  to  give  a  data  set  which  could 
be  used  to  characterize  the  explosive  test  seismically.  Simple  cube  root  of 
yield  scaling  was  used  to  give  array  dimensions  comparable  to  the  2-20  km 


spaclngs  we  have  used  previously  for  nuclear  shots  at  MTS.  A  surface  view  of  the 
near-field  array,  along  with  some  of  radial  close-in  gages.  Is  given  In  Figure 
3.  The  near-field  array  spanned  the  50  to  150-meter  rangfe  with  the  largest  azi¬ 
muthal  separation  of  gages  being  60°.  Each  station  consisted  of  a  3-component 
force  balanced  accelerometer  recorded  by  the  DR-100  digital  event  recorder. 

The  data  were  sampled  at  200  samples  per  second. with  a  5-pole  butterworth 
filter  at  50  Hz.  Eleven  stations  were  fielded  for  a  total  of  33  possible  - 
records. 

The  final  instruments  were  placed  at  1.6  km  north  and  west  of  the  shot 
point.  From  previous  tests,  it  has  been  learned  that  at  these  distances,  fun¬ 
damental  and  higher  made  surface  waves  become  important.  Each  station  con¬ 
sisted  of  a  3-component  Sprengnether  S-6000  velocity  gage  recorded  by  a  Terra 
Technology  digital  event  recorder.  The  data  were  sampled  at  200  samples  per 
second  with  5-pole  butterworth  filters  at  70  Hz.  Two  stations  60  meters  apart 
north  and  one  station  west  of  the  shot  were  fielded  for  a  total  of  9  com¬ 
ponents  of  motion. 
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DATA  SUMMARY 


All  the  data  recorded  from  this  experiment  are  reproduced  In  this  report. 

t 

Of  the  48  close-in  gages,  45  were  recorded,  two  of  which  clipped— gages  1334 
and  1337.  Figures  4  through  12  contain  the  close-in  data.  The  signal  quality 
and  signal  to  noise  ratios  are  good  for  the  data  set.  Longer  time  windows  are 
available  as  the  analog  tape  recorders  were  run  for  five  minutes  following  the 
explosion. 

The  near-field  3-component  accelerometer  data  In  the  50  to  150-meter  range 
are  illustrated  in  Figures  13  through  17.  Of  the  33  gages  fielded,  28  yielded 
satisfactory  data. 

Finally,  the  3-component  far-fleld  velocity  data  are  given  in  Figures  18 
through  20.  Nine  out  of  nine  gages  operated  satisfactorily,  yielding  good 
quality  data. 

Out  of  90  total  gages  fielded  in  this  experiment,  82  yielded  good  quality 
data  for  a  91  percent  return  rate  on  instrumentation.  The  signal  to  noise 
ratio  of  all  the  data  Is  high  (50  db  or  more)  and  so  will  warrant  careful 
analysis. 
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EXPLOSIVE 


The  spherical  charge  was  20*»“  in  diameter  and  was  made  of  TNT.  The  burn 

* 

rate  of  TNT  is  approximately  20,000  ft/sec  giving  a  detonation  time  of  42 
microseconds.  To  check  this  estimate,  a  number  of  similar  TNT  spheres  were 
instrumented  with  time  of  arrival  crystals  which  break  when  the  explosive 
front  reaches  them.  The  data  from  these  crystals  for  a  1000-pound  TNT  sphere 
are  summarized  In  Figure  21.  Assuming  it  takes  8.7  microseconds  to  break  the 
time  of  arrival  crystal  {supported  by  experimental  evidence),  the  burn  time 
for  the  sphere  was  55.8  microseconds  with  a  maximum  2-microsecond  variation 
across  the  sphere.  There  numbers  yield  a  burn  rate  of  24,492  ft/sec  and  a 
burn  time  for  our  253-pound  TNT  sphere  of  34.5  microseconds. 
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GEOLOGY  SUMMARY 


The  geology  of  the  alluvial  test  site  has  been  extensively  studied. 

Downhole  logging,  as  well  as  P  and  S  wave  refraction  surveys,  has  been  con¬ 
ducted  (Stump  and  Relnke,  1982).  A  summary  of  some  of  the  P  wave  refraction 
data  is  given  In  Figures  22  and  23.  By  using  this  data,  other  surveys,  geologl 
cal  hole  logs,  and  first  arrival  time  data  from  explosive  shots  at  the  site, 
the  average  velocity  structure  of  Figure  24  was  arrived  at.  A  simple  3-layer 
over  a  half -space  structure  appears  appropriate.  The  effect  of  these  layers 
will  certainly  depend  on  the  wavelength  of  the  energy  passing  through  the 
structure.  Much  of  the  observational  waveforms  indicate  that  at  some  ranges, 
a  simple  half-space  may  adequately  model  the  data. 


FUTURE  WORK 


We  feel  a  high  quality  data  set  giving  good  close-in  and  near-field 
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coverage  of  a  contained  explosive  source  has  been  created.  With  this  data 
set,  we  Intend  to  characterize  the  source  In  the  nonlinear  regime  including 
fairly  sophisticated  material  models  and  do  an  equivalent  elastic  represen¬ 
tation  of  the  source  from  a  seismic  point  of  view.  It  Is  hoped  that  the 
nonlinear  modeling  of  the  source  will  allow  us  to  understand  the  physics  of 
the  phenomena  occurring  very  close  to  the  explosive  and  that  this  physics  can 
be  used  to  explain  the  equivalent  elastic  source  representation  developed  from 
the  near-field  seismic  data. 

Comparison  of  the  modeling  and  source  characterization  of  this  scaled  chem¬ 
ical  explosion  to  similar  work  on  contained  nuclear  explosions  will  allow 
investigation  of  chemical  and  nuclear  explosive  differences  and  problems  In 
scaling.  We  are  particularly  interested  in  pursuing  the  Idea  of  scaling  the 
propagation  path  and  the  source  separately. 

In  the  long  term.  It  is  hoped  that  a  series  of  these  small  scaled  experi¬ 
ments  may  be  conducted  In  a  variety  of  geologic  materials.  This  work  should 
allow  us  to  determine  the  Influence  material  responses  have  on  the  seismic 
source  characterlzaton.  We  hope  these  experiments  could  Include  material 
asymmetries,  such  as  faults  and  bedding  planes. 
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Figure  1 
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OBSERVED  AND  COMPUTED  SEISMOGRAMS 
FROM  SMALL  YIELD  HIGH-EXPLOSIVE  EVENTS 


ROBERT  E.  REINKE 
BRIAN  W.  STUMP 

Air  Force  Weapons  Laboratory 


KENNETH  H.  OLSEN 
Los  Alamos  National  Laboratory 


Local  (1.6  km  range}  seismograms  have  been  obtained  from  seven  high- 
explosive  events  fired  In  alluvium  at  the  McCormick  Ranch* Test  'Site  a  few 
kilometers  to  the  south  of  Albuquerque,  New  Mexico.  The  recorded  events 
ranged  In  yield  from  253  pounds  to  16  tons.  The  explosive  configuration 
varied  from  above  ground  and  burled  single  charges  to  surface  bermed  and 
burled  distributed  charges.  This  series  of  shots  provided  a  unique  oppor¬ 
tunity  for  studying  the  effects  of  source  configuration,  yield,  and  geologic 
structure  upon  the  observed  seismograms. 

The  McCormick  Ranch  Test  Site  Is  situated  along  the  eastern  margin  of  the 
Rio  Grande  Valley  In  Central  New  Mexico.  The  site  Is  underlain  by  fairly  deep 
alluvium.  Most  of  the  shots  were  fired  In  near-surface  playa  deposits.  An 
Idealized  2-layer  over  half-space  geologic  profile  Is  shown  In  Figure  1.  This 
profile  represents  a  surface  layer  of  dry  alluvium  over  a  layer  of  wet  allu¬ 
vium  In  turn  overlaying  a  paleozoic  half-space.  The  properties  of  the  first 
and  second  layers  and  the  first  layer  thickness  were  determined  by  shallow 
refraction  surveys.  Properties  of  the  paleozoic  half-space  were  estimated 
from  the  results  of  a  nearby  shallow  seismic  reflection  survey. 

Figure  2  shows  a  schematic  diagram  of  the  charge  configuration  for  each  of 
the  seven  recorded  events.  Three  of  the  shots  were  single  charge  depth  of 
burst  (DOB)  events,  250  pounds  burled  2.5  meters,  253  pounds  burled  11.6 
meters  and  2.5  tons  fired  at  a  depth  of  9.8  meters.  A  fourth  single  charge 
event  consisted  of  1000  pounds  fired  surface  tangent  to  the  ground  surface. 

One  of  the  distributed  charge  events  was  a  0.5  ton  HEST  (HEST  Is  an  acronym 
for  High  Explosive  Simulation  Technique).  The  HEST  charge  configuration  con¬ 
sists  of  a  uniform  layer  of  explosive  spread  over  the  ground  surface  which  in 
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turn  Is  covered  with  an  earth  layer  In  order  to  contain  the  air  blast  and 
Increase  the  coupling  to  the  ground.  Another  of  the  distributed  charge  events 
was  a  16-ton  DIHEST  (DIHEST  Is  an  acronym  for  Direct- Induced  High-Explosive 
Simulation  Technique).  The  DIHEST  consisted  of  16  tons  of  explosive  burled  at 
a  depth  of  10  meters  In  a  row  of  14  drill  holes  spaced  2  meters  apart.  All 
charges  were  detonated  simultaneously.  The  third  distributed  charge  event  was 
a  combined  HEST-OIHEST  event  made  up  of  a  16-ton  burled  DIHEST  and  a  0.5  ton 
surface  HEST.  The  HEST  center  was  approximately  30  meters  distant  from  the 
vertical  axis  of  the  DIHEST.  Both  were  fired  simultaneously. 

On  all  events,  stations  were  placed  at  a  range  of  1.6  km  to  the  north  and 
west.  On  same  events,  a  2-component  linear  array  was  placed  at  the  north 
stations.  The  stations  were  made  up  of  3-component,  2-Hz  natural  frequency 
seismometers  recorded  on  digital  cassette  tape  at  a  sample  rate  of  200  samples 
per  second. 

Samples  of  the  records  from  the  seven  detonations  are  shown  In  Figures  3, 

4  and  5.  All  of  these  records  are  from  the  1.6  km  north  station.  The  first 
four  (250  lb  2.5  m  DOB,  253  lb  11.6  m  DOB,  0.5  ton  ST,  and  0.5  ton  HEST) 
records  in  each  figure  have  the  same  amplitude  scale  so  that  direct  com¬ 
parisons  may  be  made  among  these  twelve  records.  The  amplitude  scale  on  the 
fifth  record  (2.5  ton  9.8  m)  Is  five  times  as  large  as  the  first  four  (l.e., 
an  excursion  of  the  trace  on  record  5  represents  a  velocity  amplitude  5  times 
that  of  the  same  excursion  on  records  1  through  4).  Amplitude  scales  on  the 
last  two  records  (16-ton  DIHEST  and  16-ton  DIHEST  +  0.5  ton  HEST)  represent 
motions  25  times  greater  than  the  first  four  records  and  5  times  greater  than 
record  5. 

Examining  the  first  four  records  in  Figures  3,  4,  and  5  (the  low  yield 
events),  the  HEST  event  appears  to  be  the  most  efficient  at  coupling  energy 


Into  the  ground.  The  surface  tangent  burst  Is  perhaps  the  least  efficient  as 
far  as  energy  coupling  Is  concerned,  with  the  exception  of  the  prominent  air 
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coupled  Rayleigh  wave  appearing  at  about  4  seconds. 

The  11.6  m  DOB  records  show  a  very  nicely  dispersed  fundamental  mode 
Rayleigh  wave  train  with  a  rather  sharp  cutoff  at  the  end  Indicating  the  pre¬ 
sence  of  an  Airy  phase.  Almost  all  of  the  records  contain  two  groups  of  sur¬ 
face  waves.  The  first  arriving  group  is  likely  composed  mainly  of  higher  mode 
propagation  with  the  later  arriving  packet  being  associated  with  fundamental  mode 
propagation.  The  HEST  event  records  contain  a  third  prominent  high-frequency 
wave  group  arriving  between  the  lower  frequency  higher  and  fundamental  mode 
packets.  This  arrival  is  apparently  some  type  of  higher  mode  surface  wave  as 
well.  Since  this  particular  phase  is  not  distinct  on  the  records  from  any  of 
the  other  events.  Its  excitation  may  be  a  result  of  the  HEST  type  of  source 
configuration.  In  general,  the  fundamental  mode  surface  wave  frequency  con¬ 
tent  appears  to  vary  little  from  shot  to  shot.  The  main  changes  seem  to  occur 
In  the  higher  mode  phases. 

The  body  wave  frequency  content  changes  relatively  little  for  the  four  low 
yield  shots.  The  HEST  and  the  11.6  m  DOB  event  seem  to  be  somewhat  more  effi¬ 
cient  at  exciting  the  higher  frequencies  than  the  2.5  m  DOB  and  the  0.5  ton  ST 
shots.  The  larger  yield  shots,  as  would  be  expected,  excite  body  waves 
possessing  a  frequency  content  roughly  a  factor  of  two  lower  than  that  of  the 
four  small  yield  shots.  Unfortunately,  the  recorders  did  not  trigger  early 
enough  to  record  the  body  wave  portion  of  the  combined  HEST-DIHEST  event  so 
that  we  were  unable  to  determine  what  effect  the  simultaneous  firing  of  the 
HEST  had  on  the  lower  frequency  DIHEST  body  waves. 

Given  the  asymmetric  nature  of  some  of  the  sources,  we  had  hoped  to 
observe  differences  between  the  waveforms  recorded  at  the  north  and  west 
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locations.  Figure  6  shows  a  comparison  between  the  north  and  west  recordings 
of  three  of  the  shots.  As  seen,  there  are  Indeed  differences  between  the  two 
stations  for  the  DIHEST  event,  which  Is  a  strongly  asymmetric  source,  however, 
the  north  versus  west  differences  are  just  as  apparent  for  the  burled 
253-pound  spherical  charge.  In  general,  the  amplitudes  are  uniformly  higher 
on  the  west  for  all  shots  suggesting  that  the  observed  north  versus  west  dif¬ 
ferences  result  from  local  geologic  effects.  While  source- Induced  azfm; 
variations  may  be  present,  they  are  probably  obscured  by  the  local  geoli 
effects . 

Rayleigh  wave  group  velocity  dispersion  curves  were  computed,  using 
Haskell -Thomson  technique,  for  the  2-layer  over  half-space  model  shown  In 
Figure  1.  The  observed  dispersion  for  the  first  and  second  arriving  surface 
wave  phases  from  the  DIHEST  event,  determined  using  the  zero  crossing  method, 
was  then  plotted  (solid  circles)  on  the  theoretical  fundamental  and  first  two 
higher  mode  group  velocity  curves  (solid  lines)  shown  In  Figure  7.  While  the 
observed  group  velocity  values  are  In  the  same  range  overall  as  the  computed 
values,  the  observed  fundamental  mode  dispersion  Is  normal  while  the  computed 
fundamental  mode  dispersion  Is  inverse  (higher  frequencies  arrive  first). 

Some  adjustment  to  the  assumed  2-layer  over  half-space  model  Is  needed  to 
achieve  normal  dispersion  in  the  frequency  range  of  interest. 

A  preliminary  attempt  at  modeling  some  of  the  observed  waveforms  has  been 
made  using  the  reflectivity  code  (Fuchs  and  Mueller,  1971).  The  code  was 
written  by  Rainer  Kind  and  is  In  operation  on  the  CRAY-1  computer  at  the  Los 
Alamos  National  Laboratory.  The  initial  runs  were  made  using  the  2-layer  over 
half-space  geology  shown  In  Figure  1  and  symmetric  explosive  and  vertical 
point  force  source  functions  at  a  depth  of  10  meters.  The  experimental  wave¬ 
forms  we  attempted  to  match  for  the  first  run  were  those  from  the  HEST  and 


DIHEST  events.  The  frequency  of  the  source  function  In  the  reflectivity 
calculations  was  set  equal  to  that  of  the  first  P  wave  arrival  observed  in  the 
records  from  the  corresponding  events.  ' 


Figure  8  shows  a  comparison  between  two  of  the  calculated  vertical  wave¬ 
forms  and  the  corresponding  observed  records.  The  explosive  source  function 
was  used  for  these  calculations.  Considering  the  simplicity  of  both  the 
source  function  and  the  geologic  model  input,  the  comparisons  are  quite  good. 
The  observed  and  calculated  DIHEST  comparison  is  the  best.  There  are  a  few 
subtle  differences  between  observed  and  calculated  waveforms— the  calculated 
fundamental  mode  Rayleigh  wave  exhibits  inverse  dispersion,  the  observed  is 
normally  dispersed  (as  we  saw  in  the  calculated  versus  observed  dispersion 
comparison),  and  the  relative  amplitude  of  the  higher  mode  is  a  bit  low  in  the 
calculation  as  compared  to  the  recorded  waveform— however,  the  calculations 
capture  the  overall  character  of  the  observed  DIHEST  waveform. 

The  HEST  observed  and  calculated  comparison  is  not  as  good  as  the  jIHEST 
comparison  which  is  perhaps  to  be  expected  since  the  HEST  is  really  a  number 
of  surface  point  forces  spread  over  a  finite  area  rather  than  a  buried  explo¬ 
sive  source.  The  calculation  fails  to  pick  up  the  intermediate  high  frequency 
higher  mode  surface  wave  phase  and,  again,  the  observed  fundamental  mode 
dispersion  is  normal  while  the  calculated  is  inversely  dispersed  suggesting  a 
need  for  a  slight  change  in  either  the  velocities  or  thicknesses  of  the  input 
structure.  The  calculation  does  capture  the  higher  frequency  nature  of  the 
HEST  waves  quite  well.  Figures  9  and  10  show  a  set  of  computed  vertical 
seismograms  from  200  meters  to  2  km  for  these  two  calculations. 

Since  the  initial  modeling  phase,  a  series  of  calculations  using  shallow 
sources  and  slightly  different  frequencies  has  been  attempted.  Some  selected 
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examples  of  the  results  are  shown  In  Figures  11  through  18.  Figures  11  and  12 
are  from  a  calculation  using  a  1  meter  point  force  at  12.8  Hz.  Figures  13  and 
14  are  for  a  1  meter  explosive  source  at  12.8  Hz.  Figures  15  and  16  contain 
the  results  of  a  17-Hz  explosive  source  at  1  meter,  while  Figures  17  and  18 
are  waveforms  from  a  17-Hz  point  force  at  the  1  meter  depth.  The  explosive 
sources  appear  somewhat  higher  frequency  than  the  point  forces. 

If  the  time  scales  were  removed  from  some  of  these  records,  they  would  not 
be  unlike.  In  overall  appearance,  seismograms  recorded  at  regional  distances 
from  larger  events.  To  determine  if  the  recorded  seismograms  were  In  any  way 
scaled  versions  of  those  recorded  at  larger  ranges  from  higher  yield  shots,  a 
simple  wavelength  scaling  relationship  was  borrowed  from  the  seismic  model 
laboratory  (Goforth,  1976).  The  scaling  relationship  used  and  the  resulting 
full  scale  profile  are  shown  in  Figure  19.  For  the  sake  of  argument,  a  factor 
of  10  change  In  frequency  was  assumed  between  the  McCormick  Ranch  events  and  a 
150  kt  event  recorded  at  near  regional  distances.  Using  this  value  and 
obtaining  velocity  ratios  between  the  McCormick  Ranch  half-space  and  the  real 
earth  mantle,  an  earth  to  model  length  ratio  of  17.3  was  obtained.  The  pro¬ 
file  in  Figure  19  Is  the  full  scale  earth  profile  modeled  by  the  McCormick 
Ranch  events.  The  profile  Is  not  a  very  relevant  one  In  terms  of  any  real 
earth  structure,  however,  it  does  Illustrate  the  scaling  principle  Involved 
which  applied  to  other  test  sites  and  events  might  yield  more  realistic 
results. 

These  preliminary  results  indicate  that  body  and  surface  wave  responses 
are  fairly  well  modeled  using  the  reflectivity  technique  at  these  relatively 
close-in  distances.  The  different  source  configurations  result  In  some  subtle 
and  some  not  so  subtle  differences  In  the  observed  waveforms.  The  similarity 
in  character  of  these  waveforms  to  those  recorded  at  regional  distances  from 
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larger  events  suggests  that  It  may  be  possible  to  study  some  problems  of 
regional  wave  propagation  by  using  relatively  small  HE  events.  Although  the 
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recorded  events  and  geology  at  McCormick  Ranc>t  do  not  scale  to  any  real  earth 
situation  of  Interest,  other  sites  might  well  be  useful  for  small  scale  testing. 
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ABSTRACT 


The  proximity  of  the  MILL  RACE  test  site  to  the  Trinity  site  where  the 
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first  atomic  bomb  was  detonated  In  July  of  1945  offered  a  unique  opportunity 
to  compare  the  unusual  appearing  seismic  waves  recorded  on  the  Trinity  event 
with  those  excited  by  MILL  RACE.  Two  of  the  Leet  seismographs  used  to  record 
the  1945  Trinity  event  were  obtained  and  fielded  on  MILL  RACE  alongside  modem 
digital  Instrumentation  In  an  effort  to  verify  the  Trinity  seismograms.  Due 
to  failure  of  the  Leet  Instruments,  the  experiment  was  only  partially  successful. 
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1.0  INTRODUCTION 

The  MILL  RACE  Event  was  sited  In  the  eastern  margin  of  the  Jornada  del  Muerto 
Valley,  roughly  six  kilometers  south  of  Trinity  site,  wherfr  the  first  atomic  bomb 
was  detonated  In  1945  (Figure  1).  Several  measurements  were  madt  of  the  strdng 
seismic  motions  excited  by  the  Trinity  Event.  Five  Leet  3-component  strong 
motion  mechanical  seismographs  recorded  the  Trinity  motions  at  five  different 
ranges  and  azimuths  (Reference  1).  The  seismogram  from  one  of  these  stations 
(8.2  km  north  of  Trinity)  was  discussed  by  Leet  In  a  1946  paper  (Reference  2)  and 
has  been  the  subject  of  some  controversy  since  Leet  gave  the  name  "Hydrodynamic 
Wave”  to  one  unusual  appearing  section  of  the  seismogram  because  the  particle 
motion  was  prograde  and  elliptical,  resembling  that  of  a  water  wave  (Figure  2). 

In  a  1962  paper  (Reference  3),  Leet  suggested  that  the  "Hydrodynamic  Wave" 
emanated  only  from  explosions  and  might  be  useful  In  discriminating  between  the 
seismic  waves  produced  by  nuclear  explosions  and  those  resulting  from  naturally 
occurring  earthquakes. 

The  DICE  THROW  Event,  fired  a  few  kilometers  west  of  Trinity  in  1976,  was 
recorded  by  a  fairly  extensive  array  of  close-in  seismic  stations  utilizing  modern 
analog  and  digital  Instrumentation.  Some  of  the  seismograms  (Figure  6)  recorded 
from  this  array  closely  resembled  Leet's  1945  Trinity  record  (References  4  and  6). 
Reference  4  attempted  to  explain  Leet's  "Hydrodynamic  Wave"  and  the  similar 
appearing  waves  observed  on  the  DICE  THROW  Event  as  higher  mode  Rayleigh  surface 
waves.  The  availability  of  the  original  Leet  Instruw.ts  for  the  MILL  RACE  Event 
offered  a  unique  opportunity  to  once  again  record  the  "Hvdrodynamlc  Wave"  on  the 
original  Trinity  Instrumentation.  By  placing  modern  digital  Instruments  along¬ 
side  the  Leet  seismographs,  it  was  hoped  that  the  Trinity  seismograms  could  be 
verl fled. 
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2.0  DESIGN  OF  THE  EXPERIMENT 

/ 

2.1  The  Leet  Seismographs. 

The  Leet  Instruments  were  designed  and  constructed  primarily  for  the  regis¬ 
tration  of  vibrations  from  dynamite  blasts,  traffic,  machinery,  and  general  Indus¬ 
trial  sources.  The  Instruments  have  roughly  a  magnification  factor  of  25  for 
frequencies  above  3  Hz.  The  seismograph  Is  optical  and  mechanical  In  operation 
in  that  mirrors  are  attached  to  each  of  the  three  Inertia  elements.  Light  frou.  a 
single  filament  galvanometer  lamp  Is  reflected  from  the  mirror  and  onto  a  moving 
strip  of  photographic  paper  after  being  reflected  several  more  times  In  order 
to  lengthen  the  optical  path  and,  hence,  increase  the  magnification  (Figure  3, 
Reference  5). 

A  total  of  five  Leet  3-component  strong  motion  seismographs  were  set  out 
to  record  the  strong  motion  seismic  waves  resulting  from  the  Trinity  Test.  In 
addition  to  the  one  at  8.2  km  north,  they  were  also  placed  In  the  neighboring  towns 
of  San  Antonio,  Carrlzozo  and  Tularosa  and  at  Elephant  Butte  Lake  (Reference  1). 

Two  of  these  original  Instruments  were  found  to  be  In  the  archives  of  the  Los 
Alamos  National  Laboratory  (LANL).  LANL  archive  personnel  were  kind  enough  to 
grant  permission  to  field  the  Leet  Instruments  on  the  MILL  RACE  Event." 

The  Leet  Instruments  were  found  to  be  In  fairly  good  condition  after  some 
36  years  In  storage,  with  the  exception  of  some  missing  parts.  Missing  from  one 
Instrument  were  several  gears  from  the  camera  drive  train  and  the  camera  itself. 

The  other  instrument  was  complete  and  in  operating  condition.  Using  the  complete 
instrument  as  a  model,  LANL  technicians  were  able  to  fabricate  the  camera  and 
drive  train  sections  missing  from  the  second  instrument. 

The  instruments  were  then  bench  tested,  whereupon.  It  was  discovered  that 
the  film  moved  past  the  camera  slit  at  a  speed  of  only  about  two  Inches  per  minute 
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rather  than  the  average  speed  of  about  four  Inches  per  second  mentioned  by  Leet 
In  Reference  5  and  roughly  the  apparent  film  speed  of  the  Instrument  which  pro¬ 
duced  the  record  from  Trinity  (References  1  and  2).  Apparently,  the  two  archive 
Instruments  were  used  at  the  more  distant  ranges  on  Trinity  where  a  longer 
recording  time  was  needed,  and  resolution  of  the  higher  frequencies  was  not  so 
important.  To  remedy  the  problem,  the  LANL  geophysics  division  procured  two  higher 
RPM  motors  which  were  installed  by  AFWL/NTESC  technicians,  bringing  the  average 
film  speed  to  near  four  Inches  per  second. 

Several  additional  problems  were  encountered  during  the  fielding  of  the  Leet 
Instruments.  The  seismographs  were  designed  to  operate  on  110  volt,  60  cycle  A.C. 
power  which  was  not  available  at  the  desired  station  location.  This  necessitated 
the  use  of  two  24-volt  D.C.  to  110-volt  A.C.  inverters  to  provide  the  required 
power.  Since  safety  requirements  associated  with  the  MILL  RACE  Drone  experiment 
prohibited  the  manning  of  the  Leet  seismograph  station,  a  way  was  needed  to  turn 
the  Leet  instruments  on  shortly  before  shot  time.  Initially,  It  was  decided  to 
use  timers  to  apply  the  24-volt  power  supply  to  the  Inverters  approximately  30 
seconds  before  shot  time.  This  plan  was  ruled  out  when  It  was  learned  that, 
because  of  the  drone  experiments,  shot  time  would  not  be  known  In  advance  to  within 
two  or  three  minutes. 

The  lack  of  advance  knowledge  of  the  shot  time  made  it  necessary  to  find  a 
way  to  turn  on  the  Leet  Instruments  remotely  from  the  observation  point.  To 
accomplish  this,  AFWL/NTESC  modified  a  radio-controlled  firing  system  normally 
used  for  firing  explosive  charges  while  conducting  seismic  refraction  surveys. 

The  system  was  modified  so  that,  rather  than  detonating  a  charge,  24-volt  power 
was  switched  on  to  the  Inverters,  which  In-turn  supplied  110-volt  A.C.  power  to 
the  Leet  instruments. 
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2.2  Modem  Instrumentation  Used  and  Instrument  Locations'. 

Two  station  locations  were  selected  for  the  MILL  RACE  Event  (Figure  1). 
One  was  approximately  8  km  slightly  west  of  north  of  MILL  RACE;  the  other  was 
approximately  8  km  west-northwest  of  the  MILL  RACE  Test  Bed.  The  northern  loca¬ 
tion  was  approximately  along  the  same  radial  from  Trinity  as  was  the  Leet 
seismograph  station  on  the  Trinity  Event.  It  was  of  course  much  nearer  Trinity 
than  the  original  1945  station,  but  was  at  approximately  the  same  range  (8  km) 
from  MILL  RACE  as  the  original  station  was  from  Trinity. 

Both  Leet  seismographs  were  placed  at  the  north  station.  Two  modern 
digital  event  recorders  were  placed  alongside  the  Leet  Instruments  for  comparison 
purposes.  These  Instruments  record  digital  data  directly  onto  magnetic  tape 
cassettes.  They  are  equipped  with  event  triggers  so  that  recording  will  occur 
when  the  amplifiers  see  a  signal— If  the  definition  of  what  constitutes  a  signal 
Is  set  correctly  by  the  user.  A  3-component  2  Hz  natural  frequency  seismometer 
was  connected  to  one  digital  recorder.  Three  1  Hz  vertical  seismometers  supplied 
the  Input  for  the  other  three  channel  digital  recorder.  These  three  vertical 
seismometers  were  set  out  in  a  linear  array  with  100-meter  spacing— that  Is,  the 
three  seismometers  were  set  out  on  a  single  radial  In  Increasing  distance  away 
from  the  MILL  RACE  GZ.  The  purpose  of  this  array  was  to  obtain  closely  spaced 
waveforms  so  that  a  phase  velocity  dispersion  curve  could  be  obtained  for  the 
"Hydrodynamic  Wave"  and  other  portions  of  the  surface  wave  train. 

3.0  RESULTS 

Out  of  15  total  channels  (If  film  records  of  light  traces  can  be  properly 
referred  to  as  channels),  only  six  recorded  successfully.  The  remote  control 
system  did  successfully  turn  on  the  Leet  Instruments;  however,  neither  Instrument 
produced  a  seismogram  of  the  event.  The  film  In  one  Leet  camera  slipped  off  the 
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take-up  reel.  The  film  In  the  other  Leet  Instrument  did  run; ‘however,  no  traces 
were  visible  after  the  film  was  developed.  Each  camera  Is  equipped  with  a  slit 
and  spring  loaded  cover  for  the  silt  which  does  not  open  un£11  the  instrument 
trap  door  Is  closed  and  locked  (Figure  3).  Locking  the  trap  door  also  engages 
the  clutch  for  the  camera  drive.  The  camera  door  was  one  of  the  last  things 
checked  the  morning  of  the  shot.  This  was  checked  by  pulling  up  on  the  door 
handle  to  ensure  that  the  door  was  Indeed  locked.  In  this  Instrument,  there  was 
evidently  some  play  or  misalignment  In  the  system  which  allowed  the  door  over  the 
camera  silt  to  close  just  enough  to  block  the  light  beams.  This  situation  had  not 
occurred  In  pre-event  tests  of  the  Instrument,  evidently,  because  the  trap  door 
was  not  pulled  up  to  check  for  the  camera  drive  engagement. 

Only  two  of  the  three  digital  recorders  operated  successfully.  Selsmogtams 
from  these  (the  west  and  north  3-component  stations)  are  shown  In  Figures  4  and  5. 
The  recorder  connected  to  the  three  vertical  seismometers  did  not  trigger  and 
begin  recording  -until  arrival  of  the  alrblast.  The  average  value  of  all  three 
channels  Is  used  to  activate  the  record  trigger.  In  this  case,  the  three  vertical 
seismometers  were  placed  at  Increasing  distances  from  the  blast  so  that  the  sep¬ 
aration  between  the  seismometers  was.  In  general,  on  the  order  of  one-quarter  to" 
one-half  the  wavelength  of  the  Incoming  ground  motions,  resulting  In  destructive 
interference  and  explaining  why  the  recorder  did  not  trigger. 

Examination  of  the  two  3-component  records  obtained  from  the  north  and  west 
stations  reveals  that  amplitudes  are  somewhat  lower  than  those  observed  at  similar 
ranges  on  the  DICE  THROW  Event.  In  addition,  the  amplitudes  observed  at  the  west 
MILL  RACE  Station  sere  roughly  a  factor  of  two  higher  than  those  at  the  north 
station.  Neither  of  the  records  bears  a  strong  resemblance  tr.  'he  Trinity  or  OICE 
THROW  records  (Figures  2  and  6).  This  suggests  that  the  geologic  structure 
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underlying  the  travel  path  between  the  MILL  RACE  62  and  the  north  station, 
especially.  Is  shallower  (the  alluvial  thickness  is  less)  and  perhaps  less 
uniform  than  the  Trinity  or  DICE  THROW  GZ  to  recording  station  travel  paths. 

4.0  CONCLUSIONS 

This  experiment  was  only  a  limited  success.  Oue  to  freakish  equipment 
failures,  the  two  Leet  Instruments  originally  used  on  the  Trinity  Event  did  not 
produce  seismograms  of  the  MILL  RACE  Event.  Two  out  of  three  of  the  modem  Instru¬ 
ments  did  record  successfully.  The  digital  record  at  the  north  station,  however, 
bears  little  similarity  to  the  Leet  seismogram  of  the  Trinity  Event  or  to 
seismograms  of  the  DICE  THROW  Event  which  are  similar  to  the  Trinity  Seismogram. 
This  suggests  that  the  geologic  structure  associated  with  the  path  from  MILL  RACE 
to  the  north  station  Is  somewhat  different  from  that  along  the  path  from  Trinity 
to  the  Leet  seismograph  station.  If  future  large  high-explosive  shots  are  to  be 
fired  In  the  Trlnlty-MILL  RACE  area,  the  Leet  Instruments  should  again  be  fielded; 
however,  the  station  location  should  be,  probably,  placed  on  an  azimuth  10  or  20 
degrees  west  of  the  MILL  RACE  radial  In  order  to  move  Into  a  more  favorable 
geologic  structure  for  reproducing  the  "Hydrodynamic  Wave." 
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Figure  1.  Seismic  Station  Map  for  MILL  RACE 


RECORD  OF  EARTH  MOTION,  TRINITY  TEST,  JULY  16,  1945 

(  from  Loot,  1946) 


Leet  thru  component  portable  seismograph.  (1)  Knobs  coo* 
trailing  leveling  feet;  (2)  camera  handle;  (2)  knob  for  locking  inertia 
member ;  (4)  camera;  (f)  viewing  slit;  («)  light-tourae  boning;  (7) 
cylindrical  Ion;  (8)  half-til  vend  mirror;  (9)  adjustable  minor;  (10) 
camera-drive  dotch;  (11)  timing-line  shatter;  (12)  timing-line  light 
housing;  (13)  timing-tins  minor;  (14)  and  (l?)  plane  mirrors;  (IB) 
point  of  reflection;  (16)  permanent  magnet;  (If)  locking  not;  (19) 
spring;  (20)  and  (27)  inertia  members;  (21)  magnet  poles;  (22)  damping 
vane;  (23),  (24)  and  (26)  inertia  member  mirrors;  (26)  stabilising 
spring;  (28)  condensing  lens;  (29),  (30),  and  (31)  fixed'  minora. 


Figure  3.  Drawing  of  Leet  Seismograph 
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-20  SECONDS' 


MILL  RACE  LEET  STATION  (NORTH) 
RANGE  8.0Km 


Figure  4.  Seismograms  from  North  Digital  Recorder 


